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Abstract 


In  this  report  we  describe  progress  in  developing  and 
implementing  a  general  theoretical  approach  to  describing  the 
properties  of  defects  and  impurities  of  a  general  nature  in  non- 
metallic  solid  systems.  Also  we  review  progress  in  complementary 
non-destructive  experimental  tests.  These  tests  employ  and  extend 
newly  developed  techniques  in  EPR  and  NMR  imaging  in  solids. 

First  of  a  kind  EPR  imaging  experiments  have  been  performed,  and  a 
new  NMR  imaging  laboratory  has  been  constructed. 
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Introduction 


The  aim  of  this  research  is  to  characterize  point  and  extend 
defects  and  to  elucidate  their  role  in  initiating  chemical 
processes.  The  research  is  being  carried  out  in  two  separate  but 
interacting  branches.  The  theoretical  thrust  is  to  develop  a 
methodology  for  obtaining  absolute  energies  of  defects  (impurity 
ions,  charged  states  and  low  lying  excited  states)  in  solids. 

From  these  studies  a  model  for  the  initiation  of  energetic 
materials  is  being  formulated.  These  two  issues  are  reviewed  in 
the  first  two  sections  of  this  report.  The  first  section  outlines 
the  possible  role  of  charged  defects  in  initiating  energetic 
reactions.  The  second  section  reviews  progress  made  in  creating  a 
self  consistent  algorithm  for  calculating  total  energies  of 
defects  in  solids. 

The  experimental  thrust  of  this  program  is  to  develop, 
implement,  and  exploit  non-invasive  techniques  for  imaging  in 
solids.  The  two  experimental  techniques  that  are  being  developed 
are:  Electron  Paramagnetic  Resonance  (EPR)  and  Nuclear  Magnetic 

Resonance  (NMR) .  An  NMR  imaging  laboratory  has  been  constructed 
under  the  direction  of  one  of  the  Principal  Investigators  (B.  H. 
Suits).  An  outline  of  the  strategies  used  to  construct  the  NMR 

laboratory  is  included  in  the  third  section  of  this  report.  The 
last  two  sections  of  this  report  deal  with  the  recently  developed 
EPR  imaging  techniques.  First  of  a  kind  EPR  imaging  was  performed 
under  the  direction  of  one  of  the  Co-investigators  (S.  A. 
Marshall).  This  research  is  outlined  in  section  four  of  this 
report.  The  last  section  of  this  report  deals  with  EPR 
spectroscopy  of  magnetic  defects  in  solids. 
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Theory  Sub-pro iect  Progress 


BACKGROUND:  NH^CIO,, 
4  4 


Ammonium  perchlorate  is  a  white  crystalline  solid  which 
undergoes  a  reversible  crystallographic  transition  from  a  low- 
temperature  orthorhombic  structure  to  a  cubic  structure  at  240° C 
( 5 13° K ) .  [1-4]  This  polymorphic  change  is  attributable  to  the 

onset  of  free  rotation  of  the  perchlorate  anions.  [5,6]  The  low- 
temperature  bipyramidal  orthorhombic  form  has  cell  dimensions 
a.  =  9.202  A,  b.  =  5.816  A,  c.  =  7.449  A.  [7]  The  cubic  high- 
temperature  form  [8-9]  contains  four  molecules  per  unit  cell  which 
has  a  cube  edge  of  7.63  A.  The  calculated  density  [10]  of  the 

3  3 

cubic  form  is  1.76  g/cm  [5],  as  compared  with  1.95  g/cm  [11]  for 
the  measured  density  of  the  low-temperature  form.  The  phase 
transition  from  orthorhombic  to  cubic  is  endothermic  and  is 
accompanied  by  a  significant  amount  of  decomposition  and  heat 
generation.  [12]  The  computed  heat  of  transition  form  the 
orthorhombic  to  th  cubic  crystal  structure  is  2.3  ±  0.2  kcal/mole. 
[10]  There  is  some  uncertainty  as  to  the  effect  of  the  phase 

\ 

transition  on  the  kinetics  of  ammonium  perchlorate  decomposition. 
Bircumshaw  and  Newman  [13]  first  observed  that  the  maximum 
decomposition  rate  rose  sharply  to  a  maximum  at  238°,  then  fell  to 
a  minimum  at  250°,  and  finally  increased  again.  Rate  constants 
for  the  cubic  form  were  substantially  lower  than  those  for  the 
orthorhombic  form  Shidlovskii,  et  al. ,  [14]  and  Manelis  and 

Rubtsov  [15]  have  also  reported  that  the  kinetic  description  of 
the  decomposition  curves  depends  on  the  crystal  form.  In  the 
decompositon  of  the  cubic  modification  above  240°  decomposition 
centers  are  initiated  by  growth  of  a  single  nucleus  (rather  than 


»*«  V, 


by  aggregation),  are  spherical  in  form,  and  are  randomly 
distributed  throughout  the  body  of  the  crystal.  The  activation 
energy  for  nucleus  growth  is  17  ±  1  kcal/mole.  [12]  There 
aparently  has  been  little  work  done  on  the  role  of  dislocations  in 
the  decomposition  of  ammonium  perchlorate.  The  effect  of  heating 
ammonium  perchlorate  to  the  decomposition  temperature  is  to 
increase  the  dislocation  density  in  the  vicinity  of  a  nucleus. 
Least-mean-square  values  for  the  activation  energies  are  33.91  ± 
1.64  kcal/mole  for  orthorhombic  and  26.97  ±  2.36  kcal/mole  for  the 
cubic  modification.  [16]  The  question  of  the  precise  effect  of 
the  phase  transition  on  the  decomposition  is  still  an  open  one. 

Bircumshaw  and  Newman  [17,18]  considered  three  possible 
mechanisms  for  the  decomposition  of  ammonium  perchlorate.  They 
are 

1.  Electron  transfer  from  a  CLO^  anion  to  an  (interstitial) 
NH*  cation. 

4 

2.  Proton  transfer  from  NhT  to  CLO . . 

4  4 

3.  Thermal  breakdown  of  CLO^  anions  by  rupture  of  a  CL-0 
bond . 

They  ascribed  the  high-temperature  reaction  to  (3),  sublimation  to 
(2),  and  thelow-temperature  reaction  to  (1). 

Jacobs  and  Russel 1 -Jones  [19]  however  indicate  that  a  single 

unified  mechanism  is  operative  in  both  the  low-  and  high- 
temperature  decomposition  and  for  sublimation.  That  unified 

mechanism  they  believe  to  be  proton  transfer.  They  represent  the 

various  modes  of  decomposition  of  ammonium  perchlorate  by  the 

following  scheme. 


(2) 


wri.i  |.i 


NHf  <*>  +  R  -l 

NHtcLO"  <====>  NH  (a)  +  HCLO^(a) 

4  4  ?  t 


R  --->  P 


sublimate  -  NH-  (g)  +  HCLO  (g) 

III  f  T 


NH3  (g)  + 


R  fr>  P 


I  denotes  the  Low-temperature  mechanism,  II  the  high-temperature 


mechanism,  and  III  the  sublimation  process.  (a)  denotes  the 
adsorbed  state  and  (g)  the  gas  phase.  R  denotes  radical 
intermediates  and  P  the  final  products.  The  rate  of  decomposition 
is  faster  than  the  sublimation.  [12]  The  low-temperature  reaction 

must  therefore  involve  adsorbed  NH_  and  HCLO..  Early  observations 

3  4 

of  nucleation  [20]  indicate  that  certain  crystal  sites  are 
preferentially  involved.  This  indicates  the  role  of  impurities  or 
of  imperfections  in  initiating  this  reaction.  Since  the  low- 
temperature  reaction  is  suppressed  in  specially  purified  ammonium 
perchlorate,  [21]  impurities  rather  than  dislocations  [22]  are 
indicated  as  nucleus-forming  sites.  If  the  temperature  of  the 
residue  is  raised,  NHg  desorbs  from  the  crystal  surface  and  the 
high-temperature  reaction  sets  in.  The  distinction  between  this 
reaction  and  sublimation  is  whether  or  not  the  HCL04  decomposes  in 
the  gas  phase.  The  situation  gets  markedly  more  complicated  when 
catalysts  and/or  inhibitors  and  polymeric  fuels  are  involved.  In 
general  catalysts  lower  the  minimum  ignition  temperature.  A 
unified  treatment  of  the  mathematical  theories  of  ignition  is 


available  in  a  survey  article  by  price,  et  al.  [35] 


The  unified  mechanism  [19]  proposed  by  Jacobs  and  Russel  1- 
Jones  (equation  2)  has  the  NH4+CL04  going  to  neutral  NH3  and 
HCLO^ .  Isolating  this  change  of  state  and  structure  was  explored 
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by  Cardy,  et  al .  132]  They  make  a  theoretic  prediction  of  the 
kinetic  parameters  of  the  radical  formation  and  dissociation 
reactions.  The  ground  state  of  the  ammonium  radical  appears  to  be 
of  the  Rydberg  type.  They  found  its  ionization  potential  to  be 
4.29  eV. ,  the  NH^  formation  reaction  from  NHg  +  H  to  be  slightly 
endothermic.  They  predict  a  lifetime  for  NH4  of  0.1  jjs. 

Ng  and  Jacobs  in  their  paper  on  thermal  depolarization  in 
ammonium  chloride  crystals  [36]  propose  that  protons  transfer  from 
the  NH^  to  form  the  molecular  defect,  NH0  +  HCL.  The  HCL  is  then 
mobile  in  the  anion  sublattice  because  of  the  mobility  of  the 
trapped  protons.  At  the  surface  HCL  molecules  escape,  leaving 
behind  anion  vacancies,  which  diffuse  into  the  crystal  and 
eventually  trap  the  mobile  proton  holes.  Anion  vacancies  are 
known  to  be  mobile  from  diff fusion  experiments.  HCL  diffuses  down 
grain  boundaries  and  dislocations  as  HCL  molecules,  and  then  into 
the  crystal  as  protons  in  the  anion  sublattice.  Once  the  protons 
annhiliate  the  proton  holes  the  anion  vacancies  liberated  diffuse 
to  dislocations  and  anhililate  the  cation  vacancies  produced  when 
HCL  enters  the  lattice.  Perhaps  this  or  a  similar  mechanism  is 
operative  in  ammonium  perchlorate  also. 

Harris  and  Fiasson  [37]  in  their  paper  on  the  measurement  of 
surface  and  boundary  charge  on  NaCL  state,  "The  difference  between 
the  Gibbs  free  energies  required  for  anion  an  cation  vacancy 
formation  in  an  ionic  crystal  produces  and  excesss  of  vacancies  of 
one  sign  at  source  and  sinks  of  vacancies,  and  hence  causes 
dislocai tons ,  surfaces  and  grain  boundaries  to  become  charged." 
With  a  relatively  simple  apparatus  they  determined  the  potential 
on  freshly  cleaved  surfaces  of  nominally  pure,  single  crystal  NaCl 
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to  be  negative  at  room  temperature.  They  found  the  potential  near 

iS! 

» 
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grain  boundaries  in  bicrystals  to  be  positive.  They  attribute  the 

•V 

difference  to  grain  boundary  precipitation  and  predict  a  boundary 

H 

isoelectric  point  at  at  higher  temperatures.  The  presence  of 

Hi 

charges,  both  positive  and  negative,  can  be  predicted  for  the 

J 

8 

ionic  ammonium  perchlorate  crystal. 

I 

The  existence  of  the  NH^  radical,  at  least  in  the  gaseous 

r  "J 

state,  appears  no  longer  in  doubt.  The  Rydberg  Spectrum  has  been 

r  _  * 

observed  and  interpreted  by  G.  Herzberg.  [33]  He  states,  "On  the 

basis  of  the  observed  vibration  spectrum  its  generally  agreed  that 

the  NH*  ion  has  tetrahedral  structure  (point  group  T , ) .  If  an 

4  a 

electron  is  added  in  a  Rydberg  orbital  one  may  expect  the  same  T^ 

>3 

*>* 

structure  for  the  resulting  neutral  NH4  radical  in  its  Rydberg 

$ 

states . " 

% 

The  ab  initio  calculations  of  Kaspar,  Smith  and  McMaster  [34] 

on  the  potential  barrier  for  the  dissociation  of  the  ammonium 

‘-.j 

radical 

NH* - >  NH3  +  H+  (3) 

4 

h 

proved  very  helpful  to  this  work.  Their  published  values  of 

>5 

I . 

56.55903  hy .  for  uncorrelated  energy  and  -56.77444  hy .  for 

."v.l 

correlated  system  energy  of  NH^  compares  favorably  with  our  free 

radical  energies  of  -56.403662  hy.  and  -56.506768  hy. 

5? 

respectively.  Their  calculations  employed  Slater  type  orbitals 

.  .N 

and  Cl  energy  corrections.  Their  Cl  correlation  accounts  for 

.‘S 

0.221541  Hy .  or  56.86  eV.  compared  to  0.1031455  Hy .  or  2.81  eV . 

for  our  MBPT  correction.  The  difference  between  their  calculated 
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value  at  the  UHF  level  and  ours  was  0.15541  Hy.  or  4.23  eV.  The 
difference  at  the  correlated  level  was  0.26767  Hy.  or  7.28  eV. 
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THEORY  -  UHF 

Ab  initio  unrestricted  Hartree-Fock  linek 


r  combination  of 


atomic  orbital  (LCAO)  theory  is  the  method  used  in  this  study  to 
calculate  energy  values.  We  are  primarily  interested  in  the 
electronic  characteristics  of  molecular  solids.  The  problem  can 
be  divided  into  two  parts,  that  involving  motion  of  nuclei  and 
that  involving  the  electrons.  [24] 


H  t  =  E  Y 


is  the  Schroedinger  Wave  Equation.  The  hamiltonian,  H,  operates 
on  the  w function,  Y ,  giving  the  total  energy,  E.  E  is  the  sum 
of  the  potential  energy  and  kinetic  energy  of  the  system.  This 
hamiltonian  can  be  written  as 


H  =  H  .  +  H  =  H  . 

el  nuc  el-nuc 


el  represents  the  electron  contributions,  nuc  the  contributions 
from  interactions  of  the  different  nuclei  with  each  other,  and  el- 
nuc  the  electron-nuclei  portion.  Effects  from  the  external 
environment  are  neglected. 


"el  M, 


H  -  -  Z 
nuc  .  . 

i  =  l 


1  n 


1  „  1  N  Z  Z 

2  j  k 

- - 7 .  + - Z  -- - 

2M.  Vi  2  jkk  Rjk 


N  N 


el-nuc 


=  "k^l  it 


1  |  R,  -  r  .  | 

1  k  l 1 2 


n  is  the  number  of  electrons.  r  describes  the  electron 


coordinates.  N  is  the  number  of  nuclei,  R  their  coordinates,  M 


their  mass,  and  Z  their  charge.  The  nuclei  are  assumed  infinitely 
massive  and  the  Born-Oppenheimer  (adiabatic)  approximation 
employed.  The  calculations  then  involve  only  the  average 

o 

positions  of  the  nuclei  [24]  and  H  ,  is  redefined  as  H  , 

e  el-nuc  el-nuc 

The  equation  for  the  motion  of  the  electrons  then  becomes 


(H  .  +  H  .  )  v|/(r)  =  E  ,  *(r) 

el  el-nuc  el  v  ' 


This  is  the  electronic  Schroedinger  wave  equation.  It  is  assumed 
the  total  wave  function,  Y ,  can  be  approximated  by  a  product  of 
nuclear  and  electron  wavefunctions .  [25]  The  wave  function  is 

redefined  as  vp  and,  it  will  be  shown,  is  a  function  of  both 
position  and  spin. 

Hex tree-Fock  theory  approximates  the  wave  function,  as  a 
finite  linear  combination  of  Slater  determinants.  [26] 

Y  =  I  C .A  (10) 

i  1  1 

A ^  is  a  Slater  determinant  and  is  a  convenient  way  to  invoke  the 
Pauli  Exclusiton  Principle.  We  use  only  a  single  Slater 
determinant  which  consists  of  products  of  spin  orbitals,  4>  .  . 

♦  j(l)  •  •  -<t>n(l) 

4>j  (n)  •  •  .4>n(n) 

l//n!  is  for  normalization,  n  is  the  number  of  electrons.  The 
determinant  allows  for  all  possible  permutations  of  electrons, 
each  electron  being  indistinguishable  from  any  other  is  equally 
likely  to  occupy  any  available  spin  orbital.  Note,  exchanging  two 
electrons  interchanges  two  columns  resulting  in  a  sing  change  for 


(11  ) 


the  determinant.  The  wave  function  is  said  to  be  antisymmetric 
with  respect  to  the  interchange  of  two  particles.  Also,  if  two 
electrons  have  the  same  space  and  spin  coordinates,  two  of  the 
columns  are  identical,  the  determinant  equals  zero;  ie.  the 
wavefunction  vanishes. 

The  unrestricted  version  of  hartree-Fock  theory  (UHF)  has  no 
symmetry  restrictions  imposed  on  <J>  other  than  they  be  eigenstates 
of  sz .  A  single  determinant  of  doubly  occupied  orbitals 
represents  a  totally  symmetric  singlet  ground  state  of  an  atom  or 
molecule.  Apparently  most  molecular  ground  states  are  of  thsi 
type.  [23]  Singlet  ground  state  energy  changes  are  an  objective 
of  this  study. 

Theexpectation  value  of  the  energy  for  a  given  Hamiltonian  and 
wave  function  is  an  upper  limit  to  the  true  energy. 

E  =  <_e_JL  H_JL_|_>  >  e  (12, 

<  A  j  A  >  true 

A  variational  method  is  applied  to  minimize  the  energy  hence 
giving  the  best  approximation  to  the  wavefunction.  The  <t> 1  s  are 
varied  until  the  minimum  energy  is  found. 


=  0 


The  <|>j  are  constrained  to  be  orthonormal.  The  orbitals  are  not 

required  to  be  doubly  occupied  or  to  have  well  defined  symmetry. 
When  the  <J>^  are  determined  variational ly ,  they  are  determined  by 
the  Hartree-Fock  equaiton: 


F  ( p )  <*>  i  =  €  .  <t>  . 


F  is  the  Fock  operator  and  p  is  the  first  order  density  matrix. 


mm 


The  total  energy  then  becomes 


E  =  Z  e.  (15) 

i  1 

The  nuclei  are  assumed  infinitely  massive  and  the  Born-Oppenheimer 
approximation  employed. 

The  Hamiltonian  in  the  Born-Oppenheimer  approximation  is 


12  It  1  '  1 

h  -  -  -5-  e  vr  +  Z  Z  T=— ----r  +  z 

2  k  k  k  i  '  Rk  i'  2  k ,  j  I  k  rjl 


+  Z 
k,  j 


.  Z.  Z  . 

_ k__i2__ 


using  atomic  units  where  e  «  ti  =  1  and  1  hartree  =  27.2  eV.  Z. 

^  '  J 

is  the  atomic  number  of  the  nucleus  at  site  R.  ..  r.  .  .  define 

k  1  j  1 1  k  1  j 

th  electron  coordinates.  There  are  n  electrons  and  N  nuclei  in 
the  system.  The  last  term  is  a  constant  because  nuclear  motion  is 
being  neglected.  When  the  nuclei  are  assumed  fixed,  kinetic 
energy  is  zero  and  the  potential  energy  constant.  The  solution  to 
an  eigenvalue  problem  is  invariant  respecting  the  addition  of  a 
constant,  once  this  potential  energy  term  is  calculated,  it  can  be 


stored  and  added  back  into  the  total  energy  solution  when  needed. 
Equation  (8)  can  also  be  written 


H  =  Z  F  +  -  -i-  Z  9ii 

i=l  *  i*j  1J 


where  are  the  one-electron  operators  and  represents  the 

two-electron  operators,  coulomb  and  exchange.  The  energy  is  then: 


E  =  Z  <0.  (1)  |F|0 .(1)  >  + 
i  =  1 

-  [<OiU)Oj(2)|g12l<l>i(l)Oj(2)>  (18) 

-  «t>.  (l)<Dj(2)|gi2|4>j(l)<l>.{2)>) 

Use  of  the  energy  variational  principle  and  varying  E  with  respect 


to  the  <1*^,  applying  Lagraingian  multipliers  for  the  orthonormality 
constraint  on  solutions,  one  can  write  the  Hartree-Fock  equation: 


Fl4'k(1)  +  Z.  [  [  4>j(2)4>j(2)  dT 2  *k(1) 

3 

-  f  *}  (2>9l2»k«2)  dT2  V1'  ]  '  Vk'1' 


(19) 


or 


„scf  ^  _scf  ,  ,  „  _ 

H  ^k  =  Ek  ^k  20 

gcf 

where  H  is  the  self  consistent  field  Hamiltonian  operator, 
equations  are  then  solved  iteratively,  the  Hartree-Fock 
Hamiltonian  operator  being  a  function  of  the  solutions  ,<J>.  . 


The 


THEORY  -  MBPT 


Calculation  of  the  energy  corrections  is  accomplished  by 

taking  the  unrestricted  Hartree-Fock  (UHF)  solution  as  the 

starting  point  and  then  apply  the  variant  of  Raleigh-Schroedinger 

perturbation  theory,  Many  Body  Perturbation  Theory  (MBPT) .  The 

single  Slater  determinant  is  enhanced  by  perturbation  theory  to 

include  correlation  and  by  a  projection  technique  to  produce 

2 

eigenstates  of  S  . 

The  exact  Hamiltonian  can  be  represented  by  the  simplified 
Hamiltonian,  H° ,  whose  eigenvalues  and  eigenfunctions  are  known 
and  a  perturbation,  V. 


H  =  H°  +  V  so,  V  =  H  -  H° 


(21) 


The  known  eigenvalues  and  eigenfunctions  are  the  zeroth  order 
result  of  the  [JHF  self  consistent  field  calculations  in  this 
study.  They  may  be  represented  as 


,o 


ie . 


H  $  .  *  w  .  *  . 
J  J  J 


E°  =  w  .  =  <$  .1  H°|$  .> 


(22) 

(23) 


The  first  order  correction  is  the  expectation  value,  E..,  of  the 
perturbation  of  the  unperturbed  state.  Second  order  corrections 
require  that  the  wavefunction  be  correct  to  first  order. 


E1  = 
J 


<$  ,|  V|$  .>  =  V 
j  J  oo 


(24) 


addi tionally , 


The  second  order  correction  to  energy  which  we  seek  is 


(26) 


where  *  . 

J 

state  we 
unknown , 


is  the  first  order  wavef unct ion .  Defining  the  particular 
seek  as  v|>q  which  yields  energy,  Eq,  which  are  both 
then 


Applying 

write 


Hvq  =  EqYo  =  (H°  +  V)*0  (27) 

the  known  eigenvalues  and  eigenfunctions  of  H°  we  may 


(H°  -  »0)*0  »  (E0  -  V  -  »„>*„  (28) 

The  inverse  operator,  (H°-u>o)  1  ,  could  then  be  used  to  determine 

However  the  function  being  operated  upon  must  have  no 
projection  on  to  <t>Q,  otherwise  singularities  can  occur.  The 

projector,  P,  will  project  out  of  a  part  of  <J> q  so  that 
singularities  will  no  occur. 


P  =  \<t>  >«t>  | 
o  o 


(29) 


Multiplying  by  the  projection  operator,  (]  -  P) 


(1  -  P)  (E0 


(30) 


*  »,*  V.f  4. i 


(1  -  P)  (H°  -  w0)vt»0  = 


-  V  -  W0)^0 


Defining  and  applying  intermediate  normalization 


-  1 


then,  after  some  rearrangement 


=  [1  -  ( H°  -  Qq)  V  -  P)  (Eo  -  V  -  Wo)]  14»p  (32) 


This  in  turn  results  in  a  power  series  for  (j/p.  Developing  an 


expression  for  the  energy,  E 


<»0IH|»0>  =  <»0IE0I*0> 


<»0I«  l*0>  +  <0oIV|*o>  -  <»0|E0IU>0>  (34) 


Solving  yields 


E  =  w  +<<t>|V|M>> 
o  o  o '  1  o 


Replace  with  the  power  series  expansion,  construct  V  in  terms 
of  of  the  <t>j  ,  then 


V  =  <4>  I  VI  <f>  > 
oo  o'  1  o 


SSSBw 


However,  Eo  is  part  of  the  solution  and  must  be  replaced  by  its 
first  order  expression.  Let  Eq  a  w  +  Voo.  then  to  second  order 

V  .V  . 

E  =  0)  +  V  +  I  _2.2_2.1_  =  o)  +  V  +  E  (37) 

O  O  OO  .  _  <|)  -  (|)  .  o  oo  o 

J*0  o  j 

This  is  the  form  desired  for  second  order  corrections.  It  has 
been  shown  [27]  that  this  procedure,  MBPT,  gives  about  90%  of  the 
correlation  energy  available  in  a  given  basis  set. 
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EXPERIMENT 


The  ammonium  portion  of  ammonium  perchlorate  as  a  free  radical 
in  the  presence  of  a  negative  point  charge  has  been  studied  in  a 
previous  work.  [1]  The  intent  of  this  study  is  to  extend 
investigation  of  the  negative  point  charge  defect  phenomenon  to 
the  crystal  structure  and  expand  it  to  include  vacancy  type 
defects . 

Crystalline  ammonium  perchlorate  has  two  temperature  dependent 
forms.  Low  temperature  orthorhombic  and  high  temperature  cubic. 
High  temperature  cubic  was  selected  as  the  object  of  this  study. 

The  unrestricted  Hartree-Fock  calculation  utilizes  the  linear 
combination  of  atomic  orbitals  (LCAO)  approximation  to  construct 
the  spin  oroitals.  That  is 

<t> ,  =  E  c .  .  x  •  (38) 

i  aj  j 

c_  is  an  appropriately  chosen  unitary  matrix.  Xj  are  primitive 
Gaussian  type  orbitals  (GTO's). 


r  9  m  n  1  i 

=  a..  [  xA  y  z  e  J  J 


-b  .  r' 


a  .  and  b  .  are  input  parameters  for  the  POLYATOM  computer  code 
J  J 

utilized  for  calculation  of  integrals  over  atomic  orbitals.  x^ , 
ym,  and  z°  represent  the  type  of  atomic  orbital,  s,  p,  d,  f,  etc. 
Xj  are  primitive  Gaussians.  Contracted  basis  functions  used  in 
this  study  are  linear  combinations  of  primitive  gaussians.  (29] 

The  NHt  molecule  was  constructed  from  Huzinaga  atomic  basis 
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sets.  [30]  The  exception  is  hyfrogen  which  is  the  Dunning  and  Hay 
[31]  contraction  H(is)  of  the  Huzinaga  H(4s)  uncontracted  to 
H(3s/ls)  to  give  additional  freedom  to  the  system. 


All  variations  in  the  charge  array  were  made  at  the  location 


of  the  first  point  charge  on  the  negative  z-axis,  (0,0-7.0338 
a.u.),  which  was  a  negative  charge.  The  configuration  of  the 
radical  was  not  altered  from  its  tetrahedral  form  found  in  the 


orthorhombic  configuration.  The  intent  is  to  study  the 
.+ 


tetrahedral  NH4  radical  [33]  at  a  lattice  site  andthen  to  vary  the 
immediate  environment.  The  hydrogen  designated  HI  located  at 


(0,0,-1.915  a.u.)  was  moved  incrementally  down  the  z-axis.  UHF  + 


E  correlation  energy,  and  Mullikan  populations  were  monitored. 
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The  chemical  decomposition  of  a  system  such  as  NH..C10,,  is  a 

4  4 

complex  issue.  One  might  be  interested  for  instance  in  the 

removal  of  H(or  H+,H  )  from  the  NH+  unit  or  of  0(or  0  )  from  CIO. 

4  4 

and  even  in  the  formation  of  OH,  OH  ,  HgO,  HC1  etc.  In  the 
initial  case  presented  here,  we  consider  the  possible  removal  of 
1+  from  an  NH4  unit,  making  it  available  for  further  reaction. 

The  results  presented  here  are  for  the  UHF  +  MBPT  results  only. 

In  the  case  of  a  perfect  crystal.,  a  barrier  of  3.5445  eV  exists  to 
the  removal  of  a  hydrogen.  This  is  seen  in  the  figure  labled  NH* 
5x5x5  array  -  no  defects.  Here  the  z-axis  extends  from  the  N  of 
an  NH4  to  the  Cl  of  a  nearest  C104.  If  the  nearest  neighbor  Clo4 
is  removed  from  the  lattice  creating  an  anion  vacancy  (a  common 
type  lattice  defect)  the  barrier  to  H  removal  is  raised  to  5.87  eV 
as  is  seen  in  the  figure  labled  Vacancy  Defect.  In  this  case  the 
defect  has  a  net  positive  charge.  If  however  the  Cl04  is  replaced 
by  a  double  negative  ion,  yielding  a  net  negative  region  the 
barrier  to  H  removal  drops  to  1.3245  eV  (a  substantial  drop,  and 
at  distance  close  to  the  NH*  is  exothemic).  In  this  case  the  H  is 
removed  as  a  proton  this  is  seen  in  the  figure  -  chag.  defect. 

This  latter  case  is  quite  important  as  it  shows  the  presence 
of  a  net  negative  region  promoting  the  lowering  of  a  barrier  to 
chemical  change.  However  this  is  still  not  quite  realistic. 
Usually  in  a  crystal  some  from  of  overall  charge  neutrality  is 
maintained.  Thus  a  double  negative  ion  is  often  compensated  by 
the  presence  of  a  double  positive  ion  or  more  likely  by  near  by 
negative  ion  vacancy.  We  allow  for  this  as  well.  In  this  case  we 
consider  a  negative  ion  vacancy  in  one  of  three  places:  a)  in  a 


shell  thrid  neighbor  to  the  NH  ion  labled  (7,7-7)  b)  in  a  shell 
fifth  neighbor  to  the  NH*  ion  labled  (0,7-14)  c)  in  the  first 
shell  to  the  NH*  ion  labled  (0,7,0).  This  is  seen  in  th  final 
figure  labled  NH*  5x5x5  array.  This  array  is  neutral  and  consists 
of  two  types  a  perfect  crystal  which  duplicates  the  results  of  the 
first  figure  and  is  called  D=0,  and  the  double  negative  ion  with 
associated  charge  compensation.  As  is  seen  the  charge 
compensation  only  weably  modifies  the  effect  on  H+  removeal .  We 
conclude  negative  charged  regions  promote  the  availability  of  H  or 
H+  ions  for  further  chemical  reaction.  Furthermore  the  removal  of 
the  H+  leaves  a  neutial  NH^  behind  which  in  the  lattice  looks  like 
a  net  negative  region  to  a  neighboring  NH*  ion.  The  ingredients 
for  a  chain  reaction  of  sorts  may  be  present.  Currently  we  are 
completing  like  studies  for  the  C104  ion. 


Molecular  Crystal-ICECAP  (MC-ICECAP) 

The  general  procedure  called  ICECAP  has  been  described  fully 
at  a  prior  time.  The  procedure  which  is  an  Acronym  for  I(onic) 
C(rystal)  E(lectronic)  C (alculat ion )  A(automated)  P(rogram)  has 
been  given  a  reasonably  general  recent  implementation.  This 
implementation  is  tuned  to  systems  which  are  bulk  in  nature  and 
which  have  relatively  simple  unit  cells.  This  tuning  is  iptimal 
for  most  ionic  type  crystals  including  ceramics,  but  is  rather 

less  than  optimal  for  a  Van  der  Waals  type  of  system,  such  as 
nitro-methane  or  FDX  for  that  matter.  The  flaw  is  in  the 
particular  ICECAP  implementation  would  the  process  itself. 


Basically,  the  given  implementation  would  work  reasonably  well  for 


M** 


complex  systems,  in  terms  of  precision,  but  would  be  hopelessly 
slow  in  the  general  application. 


To  rectify  these  problems,  a  new  form  of  ICECAP  which  we  call 
MC-ICECAP  has  been  formulated.  The  nomenclature  refers  to  the 
essential  modif icaitons  to  the  basic  ICECAP  procedure  the  tune  it 
to  the  essential  features  of  a  Van  der  Waals  bonded  system.  There 
are  two  basic  assumptions  made  here.  The  first  is  that  there  is 
some  finite  distance  away  from  some  possible  lattice  flaw,  for 
which  the  molecules  comprising  the  system  look  the  same  as  the 
molecules  in  the  pure,  perfect,  periodic  lattice.  They  may  not  be 
sited  at  the  perfect  lattice  site,  that  is  fine,  and  they  may  be 
polarized  a  small  amount,  as  long  as  the  polarization  is  described 
by  linear  response  theory.  These  molecules  are  said  to  comprise 
the  environment.  The  remaining  molecules  termed  the  culster,  need 
not  be  governed  by  linear  response  theory.  The  second  assumption, 
and  the  one  not  made  in  normal  ICECAP  is  that  the  molecules  in  a 
Van  der  Waals  system  overlap  only  a  very  small  amount.  The  key 
ingredient  in  this  assumption  is  that  the  local  orbital  screened 
potential  is  given  in  terms  of  inter-molecular  coulomb  potentials 
alone.  These  may  be  of  a  polarization  type. 

The  essential  virtue  of  the  approach  about  to  be  described  is 
that,  the  complexity  of  a  problem  viewed  in  the  usual  band  theory 
(or  molecular  orbital)  way  is  significantly  reduced.  Consider  a 
normal  MO  picture  of  a  molecule  consisting  of  N  aotms,  each  with  n 
electrons.  If  one  uses  Hartree-Fock  theory  and  second  order 

perturbation  theory,  then  the  time  to  do  the  calculation  scales  as 
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(Nn)  ,  and  the  disc  storage  scales  as  (Nn)  .  If  the  local  orbital 
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translation  applies  and  is  made,  the  time  factor  reduces  to  N(n)  , 


i.  »*»_«  »  >•» 


and  the  disc  storage  scales  as  (Nn)*  (Nn)  .  These  potential 
savings  are  significant.  For  example  in  a  solid  such  as 
which  has  two  N202  un*ts  Per  unit  cell,  one  need  at  most  consider 
two  separate  but  connected  N202  units  independently,  rather  than  a 
cell  of  2(N202)  units. 

As  an  example,  consider  a  molecular  solid.  Assume  that  the 
molecular  orbitals  for  the  pure,  perfect,  periodic  solid  are 
known.  Consider  initially  that  the  system  consists  of  an 
environment  and  a  cluster  of  a  single  molecule.  This  restriction 
will  be  removed  later.  The  potential  due  to  the  environment  is 
given  in  terms  of  the  potentials  of  the  ions  at  each  site  N  as: 

V(R  )  =  ?  £  a")  r"i_1  +  Z  l“  Y1"  r~2  (1) 

N  t=o  m— *  *  *  m— 1  1  1 

The  coeficients  A  refer  to  the  intrinsic  potential  of  the 

environmental  molecule,  whereas  the  coeficients  I  refer  to  the 

induced  dipole  moments.  Thus  each  molecule  is  represented  as  a 

simple  1-body  potential  of  a  complicated  form.  The  other 

mqlecule,  called  the  cluster  here  is  said  to  be  a  site  a.  The 

equation  satisfied  by  the  HF  system  for  the  molecules  at  site  a  is 

a  simply: 

[F  +  l  V(R  )]  <J>  =  €<t>  (2) 

nx  A 

This  equation  which  is  often  used  by  cluster  modellers  is 
incomplete,  and  neglects  the  orthogonal ization  terms  due  to  other 

electronic  orbitals.  This  has  been  known  since  the  pioneering 

work  of  Adams  and  Gilbert  in  the  early  1960's.  More 

appropriately,  one  needs  the  Adams-Gi lbert  equation  to  describe 

the  embedded  cluster.  This  is  in  its  full  generality: 


1  .*•».« 


(3) 


[Fa  +  E  V(Rn)  =  pwp]  $  =  £  $  (3) 

However  Kun2  has  shown  (1969)  that  for  systems  with  small  overlap 
this  reduces  to  a  simpler  equation.  This  is  the  Adams-Gi lbert- 
Kunz  equation: 

[Fa  +  z  V(Rn)  -  pA  E  V(Rn)  Pa]  $  =  €  $  (4) 

NxA  N*A 

It  is  this  equation  which  we  will  use  for  MC-ICECAP. 

To  proceed  consider  a  normal  quantum  chemical  solution  method. 
Expand  the  solution  to  the  A-G-K  equation  in  a  basis  of  X^'s. 

First  one  sees 

Aii  =  <Xi|FA  1  V(V  '  ^A  L  V(V  PAlXi>  (5) 

13  1  A  N*A  N  A  N*A  N  A  j 

and 

A..  =  <x  .  I  F.  I  x  .>  +  <x.|  Z  V  (  R„)  I  x  .  > 
aj  l1  A1  j  i1  N^A  N;  j 

-  Z  Z  <x  |$  4  |  E  V{R  )4  >4  |x  > 

i  m  1  *  x  N*A  m  3 

But,  one  extracts  the  ordinary  Fock  matrix  here.  That  is 

F.  .  =  <x . I F J  x .> 
lj  A1  j 

ar\d  also  define  the  external  potential  matrix: 

V.  .  =  <x .  I  Z  V ( R  ) I x  .> 

U  N«A  N  J 

This  yeilds 

+  VU  -  ZxZn  „*A  <6> 

Define:  $  =  Z  Cnx  ;  S  =  <x  I  x  > 

n  ooo  pq  p 1  q 


<x . 4 . >  =  E  C1  S . 
i  i  p  p  ip 


So  that  finally: 

A..  =  F..  +  V.  .  -  E  E  (EC1S.  )  (ECmS  .)  (ECIECmV  )  (7) 

ij  ij  ^  ^  vp  p  ap'  vtjq  qj'  7  r2s  s  rs'  v  ' 

It  is  this  equation  that  one  solves. 

The  extension  of  this  to  more  complex  clusters  is  trivial . 

One  has  a  separate  A-G-K  equation  or  each  molecule  unit  in  the 

cluster  (hence  the  scaling  in  tine  and  storage  as  N) .  The 
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external  potential  seen  by  each  clustermolecule  includes  terras  now 
from  the  other  N-l  cluster  molecules  as  well  as  the  environment, 
but  otherwise  the  equations  solved  remain  unchanged. 

As  an  illustration  to  this  set  of  equations  consider  a  simple 
lattice  vacancy.  The  molecules  neighboring  the  vacancy  in  the 
first  shell  of  neighbors  may  well  be  alike  to  each  other  but 
different  from  the  ones  in  the  second  shell,  which  are  alike  to 
each  other,  and  in  turn  all  of  these  are  not  exactly  the  same  as 
the  molecules  at  larger  distance,  which  resemble  the  pure, 
perfect,  periodic  lattice  molecules.  Thus  for  instance  in  a 
simple  fee  solid  (such  as  Kr)  one  has  a  vacancy  at  (000),  the 
first  shell  of  surrounding  atoms  at  (110)  etc,  and  the  second 
shell  at  (200)  etc.  The  bulk  like  atoms  begin  at  (211)  etc. 
shell.  There  is  no  aotm  at  (000),  12  at  (110),  6  at  (200),  and 
infinite  number  beyond.  Neverthless  the  essential  physics 
requires  here  that  we  understand  only  3  atom.  The  coupled  A-G-K 
equations  yield  these  answers  as  advertised.  The  remaining 
solutions  are  obtained  by  use  of  symmetry.  This  is  being  coded 
for  large  complex  molecules  now.  The  computer  of  choice  for  this 
is  the  FPS  T-40  supercomputer  available  in  our  laboratory.  This 
machine  has  recently  produced  speeds  in  excess  of  273  million 
floating  point  operations  per  second  in  oru  tests,  and  is  ideal  in 
speed  and  architecture  for  these  problems.  However,  even  with  the 
speed  available,  without  the  A-G-K  process  here  described  the 
study  of  a  defect  complex  in  RDX  would  be  impossible  in  our  view. 
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NMR  Imaging  Lab 

I .  Summary 

Since  the  installation  of  the  85kG  superconducting  magnet  in 
the  summer  of  1986  a  complete  NMR  imaging  spectrometer  has  been 
contructed  and  used  to  produce  research  quality  NMR  spectra  and 
images.  Preliminary  tests  have  been  made  on  a  high  power 
sinusoidal  gradient  imaging  system  with  some  positive  results. 
Spectrometer  development  is  performed  jointly  under  this  contract 
and  under  ONR  contract  N00014-86-C-0724  subcontract  SK80D7820F. 

Preliminary  measurements  have  been  made  using  NaCl03  which 
undergoes  a  thermal  decomposition  ( 1 ) .  Some  measurements  have 
been  made  on  the  feasibilty  of  creating  defects  mechanically  in 
alkali  halides. 

II.  Details 

A.  Spectrometer  Development 
1 .  Hardware 

The  spectrometer  is  constructed  using  commercially 
available  building  blocks  interfaced  to  a  small  computer. 

The  equipment  constructed  locally  from  scratch  Includes  IEEE- 
488  ( GPIB )  interface  circuits,  a  pulse  sequence  generator,  an 
oscilloscope  display  unit,  and  a  digital  input  unit.  The 
latter  is  the  most  recent  addition  and  will  be  principally 

used  to  facilitate  interactive  processing  of  the  raw  NMR 
data.  Three  NMR  imaging  probes  have  been  constructed.  The 
third  uses  the  main  body  of  the  first  but  has  a  new  sample 

mounting  arrangement  and  a  new  set  of  gradient  coils.  At 
present  both  existing  probes  are  capable  of  static  gradient 
techniques.  One  of  them  is  also  capable  of  dynamic  gradient 
techniques  though  problems  with  the  spectrometer  itself  have 
made  attempts  at  such  measurements  unsuccessful .  Ms  are 


currently  tracking  down  the  problems  in  the  dynamic  gradient 
equipment . 

New  construction  to  install  conditioned,  isolated 
ground  recepticles  and  a  ground  bus  is  complete.  This 
improvement  in  the  electrical  system  will  aid  in  data 
collection  by  reducing  interference  from  signals  on  the 
building  ground  associated  with  normal  building  operations. 

A  new  Intermediate  frequency  system  will  be  installed 
in  the  near  future  to  facilitate  some  line  narrowing  NMR 
pulse  sequences  which  may  lead  to  increased  image  resolution. 

Software 

The  sofware  consists  of  five  basic  parts:  control  of 
the  pulse  sequence  generator  and  hence  the  experiment  to  be 
performed,  the  acquisition  of  data  with  signal  averaging, 
processing  of  accumulated  data,  the  creation  of  an  NMR  image 
from  a  set  of  data,  and  modest  processing  of  images. 

Most  of  the  software  is  written  is  fortran  with 
interfacing  and  time  critical  portions  in  assembly  language. 

Proofing  the  Spectrometer 

6  7 

NMR  spectra  have  been  obtained  for  Li  and  Li  in  an 
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aqueous  solution  and  in  an  AlCuLi  alloy,  Na  spectra  have 
been  obtained  in  an  aqueous  NaCl  solution  and  NaCl  and  NaCIO 

*3 

solid  samples.  Metallic  aluminum  filings  were  measured  using 
27 

the  A1  resonance  both  before  and  after  annealing.  NMR 

images  have  been  obtained  using  the  static  gradient  technique 
2  3 

and  the  Na  signal  from  both  powdered  NaCl  inhomogeneously 
distributed  on  a  microscope  slide  and  single  crystals  of 
NaClOg.  In  the  latter  case  the  NMR  line '•is  split  by  internal 
electric  quadrupole  fields  and  the  image  is  made  using  one  of 
the  satellites.  This  will  result  in  an  image  which  is  quite 
sensitive  to  local  lattice  distortions.  Measurements  using 
the  sinusoidal  gradient  technique  using  salt  water  achieved  a 
gradient  field  strength  of  roughly  50  gauss/cra  without 
significant  difficulties  (this  will  be  improved  with  new 
coils  under  development).  The  use  of  the  sinusoidal  gradient 


technique  to  obtain  research  quality  images  has  been  thwarted 
by  problems  within  the  spectrometer.  We  suspect  that  it  is  a 
triggering  problem  in  the  slnewave  generator  circuitry. 

Specific  Results 
General 

We  have  undertaken  the  study  of  defect  formation  in 
NaCl,  KBr,  and  NaCIO  .  The  defects  are  formed  in  the  alkali 
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halides  mechanically  whereas  we  use  the  defects  created 
during  thermal  decomposition  in  the  case  of  NaClOg.  For  NaCl 
and  NaClOg  we  can  easily  observe  the  sodium  NMR  signal, 
whereas  in  KBr  we  monitor  the  bromine  signal.  Signals  can  be 
obtained  from  chlorine  and  potassium  though  they  will  be 
weak. 


Sodium  Chlorate 

Samples  of  single  crystals  of  NaClC>3  have  been  grown 
from  an  aqueous  solution.  These  crystals  will  have  an  NaCl 
type  structure  with  some  distortion  to  accommodate  the 
chlorate  group  (1,2).  Upon  heating  above  100°C  the  samples 
lose  mass  as  the  material  begins  to  decompose.  The 
decompost ion  follows  the  general  form 


4NaC103  - >  4NaCl  +  602. 


It  is  unclear  at  present  exactly  how  the  intermediate 
steps  proceed.  It  is  hoped  that  theoretical  work  under  this 
program  will  help  to  understand  the  process  in  more  detail. 
There  is  evidence  in  the  literature  to  suggest  that  the 
decomposition  occurs  first  in  regions  of  high  stress  (1).  We 
have  measured  for  one  sample,  ground  to  a  coarse  powder, 
oxygen  evolution  using  a  mass  spectrometer.  These  results 
are  shown  in  figure  1.  Initially  air  and  hydrocarbon 
contaminants  (fingerprints,  vacuum  pump  oil)  are  seen. 

However  as  time  progresses  and  the  temperature  increases,  a 
steady  evolution  of  oxygen  is  observed  which  appears  to 
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emerge  in  bursts.  This  type  of  behavior  is  very  reminiscent 
of  the  f racto-emission  work  of  Dickenson,  et  al  (3).  After 
the  initial  few  minutes  no  sodium  or  chlorine  evolution  is 
observed.  We  have  also  measured  the  optical  absorption 
( 190nm  to  820nm  wavelength)  of  a  sample  for  possible 
comparison  with  future  theoretical  work.  The  absorption 
shows  an  exitonic  peak  typical  of  simple  solids  with  a 
maximum  at  a  wavelength  of  210  nm  (5.9  eV) . 

Imaging  experiments  were  performed  using  NaCIO  before 
and  after  heating.  The  images  which  resulted  were 
highly  distorted  and  uninterpretable.  Some  examination 
showed  that  this  was  predominantly  due  to  changes  in  the  NMR 
spectrum  with  sample  heating  by  the  gradient  coils.  We  have 
examined  this  behavior  in  a  non-imaging  mode.  Two  of  the 
resulting  spectra  are  shown  in  figure  2  which  illustrate  the 
changes  in  the  spectra  with  temperature.  These  spectra  were 
obtained  with  the  magnetic  field  oriented  near  a  cubic  axis. 
There  are  two  major  changes  with  temperature,  the  overall 
splitting  due  to  the  internal  electric  quadrupole  field  is 
seen  to  decrease  with  temperature  and  the  outer  ("satellite") 
peaks  become  broader.  These  changes  are  reversible  and  so 
are  not  directly  associated  with  defect  formation  but  a 
careful  study  may  help  to  understand  the  underlying 
microscopic  changes  in  the  unit  cell  and  hence  provide  clues 
as  to  the  decomposition  pathway.  We  plan  to  investigate  this 
further  in  the  near  future.  One  principle  feature  will  be  to 
investigate  whether  the  line  broadening  is  the  same  for  all 
nuclei  (homogeneous)  or  represents  a  distribution  of  effects 
over  the  sample  (inhomogeneous)  perhaps  due  to  differing 
degrees  of  stress  or  strain  within  the  sample.  We  will  carry 
these  experiments  to  significantly  higher  temperatures  under 
more  controlled  conditions.  Only  one  previous  Na  NMR  study 
has  been  found  which  measures  temperature  effects  (4)  and 
then  to  only  63* C. 

New  circuitry  has  been  added  to  the  spectrometer  to 
turn  off  the  gradient  coils  between  measurements  and  hence 


allow  large  gradients  with  little  heat  generation. 
Preliminary  images  show  substantial  improvements  (fig.  3). 
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NaCl  and  KBr 

Commercially  available  crystals  of  NaCl  and  KBr  have 
been  acquired.  Preliminary  tests  with  a  drop  hammer 
constructed  locally  show  that  initially  clear  samples  can  be 
made  cloudy  when  hit  with  the  hammer.  Some  experimentation 
was  done  to  determine  the  correct  height  for  the  hammer  to 
avoid  breaking  the  samples.  It  is  likely  that  the  cloudiness 
is  due  to  microcracks  formed  in  the  sample.  As  of  this  date 
we  have  not  performed  any  conclusive  NMR  measurements  using 
these  samples. 
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Figure  1  -  Sinnary  of  mass  spectrum  analysis  of  decomposition  of 
Sodiun  Chlorate. 
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Figure  2  -  NMR  spectra  of  NaClO,  showing  the  change  in  the 
quadrupole  splitting  and  the  broadening  of  the 
qoadrupole  satellites  as  the  temperature  is  raised 
The  rero  of  the  frequency  shift  corresponds  to 
95. 18Hi*. 
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Report  on  ESR  Laboratory  Work  Performed  for  the  oNR  Contract. 
(S.  A.  Marshall) 


During  the  course  of  this  reporting  period,  a  number  of 
problems  were  initiated,  all  of  which  dealt  with  the  nature  of 
defects  in  solids.  These  problems  were  selected  partly  because  of 
their  basic  scientific  interest  and  partly  because  of  their 
potential  in  having  their  magnetic  species  used  as  labels  for 
magnetic  resonance  imaging. 

Before  proceeding  however  it  should  be  mentioned  that  some 
trials  were  made  on  systems  initially  considered  to  be  appropriate 
to  the  general  area  of  energietic  materials.  Of  these,  one  was 
determined  to  be  appropriate  to  the  long  term  goals  of  the  project 
in  that  it  experimentally  established  the  notion  of  stress 
generated  paramagnetic  species  in  an  otherwise  non-magnetic 
medium.  The  system  in  question  was  high  purity  fuzed  quartz  does 
not  exhibit  a  significant  ESR  spectrum.  That  is,  the 
concentration  of  free  radical  species  or  broken  bonds  seems 
minimal.  Subsequent  gamma-irradiation  of  this  medium  produces 
strong  ESR  signals  suggesting  that  broken  bonds  are  produced  by 
the  high  energy  photons.  This  of  course  is  by  no  means  surprising 
in  view  of  the  work  done  over  the  past  three  decades  on  radiation 
produced  defect  species  in  condensed  media.  What  is  observed  is 
the  following.  When  annealed,  high  purity  quartz  is  crushed  into 
a  fine  powder,  an  ESR  spectrum  is  observed  which  is  otherwise  not 
observed  and  which  is  in  many  respects  similar  to  the  spectrum 
observed  in  the  gamma-irradiated  medium.  This  is  believed  to  be 
an  effect  which  is  not  isolated  to  quartz  or  to  quartz-like  media 
but  may  very  well  be  one  which  takes  place  in  various  solids.  The 
reason  it  is  so  readily  observed  in  high  purity  quartz  and  not  so 
in  other  media  may  be  due  to  (a)  the  purity  of  the  medium,  (b)  its 
|  thermal  stability,  and  (c)  the  narrowness  of  its  ESR  spectral 

!  component  lines  which  are  of  the  order  of  0.5  Oersteds  which  in 

!  turn  are  so  because  of  the  fact  that  fuzed  quartz  is  a  solid  whose 

i 

atoms  are  composed  of  the  nuclei,  oxygen  (16)  and  silicon  (28), 
which  have  zero  magnetic  moments.  It  is  suggested  that  work  along 
these  lines  be  continued  in  view  of  its  potential  relevance  to  the 

t 
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study  of  instabilities  of  energetic  materials  to  mechancial 
stress . 

In  another  area,  that  of  the  response  of  single  crystals  of 
magnesium  oxide  (periclase)  to  mechanical  stress,  a  study  due  to 
J.  T.  Dickinson,  L.  C.  Jensen,  and  M.  R.  Mckay  has  shown  that  in 
some  specific  "cloudy"  specimens  of  magnesium  oxide,  the  fracture 
induced  emission  of  gases  was  found  to  be  significantly  higher 
than  in  "clear"  specimens  of  this  material.  In  view  of  these 
results,  we  have  undertaken  an  ESR  study  of  both  "clear"  (Horton 
Company)  MgO  and  "cloudy"  (Super  Refractories)  MgO .  The 
abbreviated  results  and  conclusions  of  our  study  on  these  two 
specimens  of  magnesium  oxide  are  a  result  of  evaluations  of  the 
spectrum  of  the  divalent  manganese  ion  impurities.  Both  crystals 
exhibit  the  spectrum  of  divalent  manganese  which  nominally 
consists  of  thirty  allowed  lines.  Furthermore,  we  find  for  both 
crystals,  that  the  manganese  ESR  spectra  yield  identical  data  on 
the  ion's  magnetic  moment,  its  nuclear  hyperfine  interaction 
energy  as  well  as  on  it  s  cubic  crystalline  electric  field  energy. 
However,  the  two  crvstals  demonstrate  some  unusual  differences. 

For  one,  the  M  =  +1/2  to  +3/2  transitions,  although  observable, 
were  much  broader  than  those  in  the  "clear"  crystals.  Our 
interpretation  of  these  observations  is  that  in  both  the  "cloudy" 
and  "clear"  crystals,  a  mosaic  effect  is  operating  and  that  in  the 
"cloudy"  specimen  it  far  exceeds  the  effect  which  takes  place  in 
the  "clear"  specimen.  From  our  preliminary  studies  we  have  been 
able  to  identify  as  well  as  quantified  this  effect.  We  find  that 
in  the  "cloudy"  MgO  specimen  there  exists  a  mid-range  disordering 
that  takes  the  form  of  a  mosaic  structure  whose  Gaussian  width  is 
5.0  degrees-of -arc  and  whose  crystalline  electric  field  amplitude 
is  uncertain  by  aabout  1.0  per  cent  while  in  the  "clear"  specimen, 
the  mosaic  width  is  under  0.50  degree-of -arc  and  its  crystalline 
electric  field  uncertainty  is  of  the  order  of  0.1  per  cent.  We 
therefore  conclude  that  the  two  crystals  are  basically  the  same 
except  that  in  one,  the  "cloudy"  crystal,  there  exists  a  disorder 
which  has  a  significant  effect  upon  spectral  line  wideths  but  on 
average  leaves  spectral  components  unaffected.  In  a  effort  to 
understand  the  problem  a  bit  more  fully,  an  electron  microscopy 
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study  of  these  two  crystals  was  undertaken  by  an  under  graduate 
student  of  the  MTU 1 s  physics  department.  This  study  proved  to  be 
very  informative  in  that  it  demonstrated  a  possible  source  of  this 
disorder.  This  microscopy  study,  although  as  yet  incomplete, 
demonstrates  the  existence  of  voids  in  the  "cloudy"  specimen.  (In 
fact,  no  evidence  of  voids  in  the  "clear"  magnesium  oxide  could  be 
detected  from  the  electron  microscopy  photographs.)  These  voids 
are  identified  as  gaseous  inclusions  trapped  perhaps  during  the 
arc-fusion  growth  process. 

A  third  study  undertaken  during  the  course  of  this  year's 
effort  was  that  of  the  superhyper fine  structure  spectrum  of  the 
divalent  manganese  ion  in  single  crystal  zinc  selenide  having  the 
zinc  blende  phase.  The  object  of  this  study  was  to  determine  how 
nearly  perfect  a  single  crystal  could  be  obtained  of  a  II-VI 
compound  and  to  then  study  the  ESR  spectrum  of  a  paramagnetic 
impurity,  say  divalent  manganese,  before  and  after  stresses  of 
various  kinds  are  applied.  We  have  completed  the  first  phase  of 
this  work  and  have  measured  with  some  precision  the  following:  (a) 
the  self  hyperfine  structure  energy  of  manganese  (55),  (b)  the 
first  shell  superhyperf ine  structure  (SHFS)  energies  of  the 
selenium  (77)  ion,  (c)  the  second  shell  SHFS  energies  of  zinc 
(67),  and  finally  the  third  shell  SHFS  energies  of  selenium  (77). 
This  is  a  rich  spectrum  and  should  be  capable  of  being  used  to 
provide  information  of  the  nature  and  degree  ofdisorder  introduced 
by  various  stresses. 

During  the  course  of  this  one  year's  effort,  two  oral  papers 
were  presented,  one  at  the  Rocky  Mountain  Magnetic  Resonance 
Symposium  held  in  Denver,  Colorado  on  Imaging  of  Electromagnetic 
Fields  in  Microwave  Resonators,  and  one  at  the  1987  Midwest  Solid 
State  Conference  held  in  Notre  Dame  University,  Norte  Dame, 

Indiana  on  The  Superhyperf ine  STructure  Spectrum  of  Divalent 
manganese  In  Single  Crystal  Zinc  Selenide.  In  addition  two  papers 

were  submitted  and  accepted  for  publication  in  1988,  one  on 
Angular  Momentum  Operator  Equivalents  and  one  on  Imaging  of 
Electromagnetic  Fields  in  Microwave  Resonators.  Further,  one 
paper  has  been  submitted  for  publication  on  The  ESR  Spectrum  of 
Tetravalent  Chromium  in  Single  Crystal  Al f a-Alumina .  Also,  two 


papers  which  bear  upon  the  work  performed  on  this  contract  are 
currently  being  drafted  for  future  publication.  These  deal  with 
the  ESR  spectra  of  manganese  in  "cloudy"  and  "Clear"  magnesium 
oxide  and  another  on  manganese  in  zinc  selenide. 

Finally,  one  graduate-  student,  Mr.  M.  T.  Umlor,  received  the 
Master  of  Science  degree  in  physics  and  wrote  a  thesis  on  work 
performed  during  the  course  of  this  contract.  it  should  also  be 
noted  that  the  following  people  made  contributions  to  this 
research  effort:  G.  Dace,  and  MTU  undergraduate,  C.  Lai, 

Y.  N.  Zhang,  M.  T.  Umlor,  all  MTU  graduate  students,  Professor  D 
Yoder-Short  of  MTU  and  Professor  J.  Furdyna  of  Nortre  Dame 
University . 


A.  Barry  Kunz 
Professional  History 

Dr.  A.  Barry  Kunz  was  born  in  Philadelphia,  Pennsylvania, 
October  2,  1940,  to  Commander  and  Mrs.  Albert  Kunz.  He  received 
his  early  education  in  the  Allentown,  Pennsylvania  Public  School 
System.  He  received  his  B.S.  degree  cum  laude  in  physics  from 
Muhlenberg  College  in  1962.  He  entered  Lehigh  University  in  the 
fall  of  1962  to  pursue  graduate  studies  in  physics.  He  received 
his  M.S.  Degree  in  January  1964  and  his  Ph.D.  in  June  1966.  His 
graduate  studies  were  aided  by  a  graduate  fellowship  from  NASA. 
From  1966  to  1969  he  remained  at  Lehigh  University  as  an 
instructor  in  physics  and  research  associate  on  projects  supported 
by  the  Atomic  Energy  Commission  and  the  Air  Force  Office  of 
Scientific  Research.  In  September  1969  he  became  a  Research 
Assistant  Professor  of  Physics  at  the  University  of  Illinois  at 
Urbana-Champaign  and  a  member  of  the  Materials  Research 
Laboratory.  He  received  several  promotions  including  that  to 
Professor  of  Physics  in  1976,  a  post  he  held  until  August  1984. 
During  the  period  1971-1974  he  was  employed  part-time  by  the  Solid 
State  Section  of  the  Aerospace  Research  Laboratory,  WPAFB,  in 
Dayton,  Ohio  as  a  research  analyst.  In  1973  he  was  co-director  of 
the  NATO  Advanced  Study  Institute  on  "Elementary  Excitations"  in 
Antwerp,  Belgium.  He  was  a  consultant  to  the  photoproducts 

division  of  E.  J.  Dupont,  Inc.  from  1974-1979.  He  has  had  an 
appointment  to  the  AFOSR,  physics  research  evaluation  group. 

During  the  calendar  year  1981  he  was  visiting  professor  at  MTU. 
Currently  he  is  Presidential  Professor  of  Physics  at  Michigan 
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Technological  University  and  head  of  the  Physics  Department.  In 
1981  he  organized  the  "Frontiers  of  Physics  Symposium"  at  MTU, 
which  included  Nobel  Laureate,  Willis  Lamb  and  Buckley  Laureate, 
Albert  Overhauser  among  the  speakers.  He  has  served  on  several 
organizing  committees  for  the  Sanibel  Symposia  organized  by  the 
University  of  Florida.  He  maintains  working  relationships  with 
several  organizations  external  to  UIUC  or  MTU.  These  include  S. 

B.  Trickey  and  R.  Bartlett  at  the  University  of  Florida,  J.  M. 
Kauffman  at  the  Johns  Hopkins  University,  and  C.  P.  Flynn  at  the 
University  of  Illinois. 

He  is  a  consultant  to  Los  Alamos  National  Laboratory  and  to 
Energetics,  Inc.  In  1986  he  was  awarded  the  MTU  distinguished 
research  award  and  in  1987  the  Michigan  Universities  Board  of 
Governor's  outstanding  faculty  award.  He  serves  on  the  editorial 
board  of  the  International  Journal  of  Quantum  Chemistry  and  of 
Computers  in  Physics.  He  is  a  Fellow  of  the  American  Physical 
Society,  a  member  of  Sigma  Xi ,  Phi  Alpha  Theta  and  Sigma  Pi  Sigma 
honorary  societies. 

He  has  been  a  pioneer  in  the  implementation  of  Hartree-Fock 
crystal  code  and  results  (LiCl).  In  1971  he  extended  the 
electronic  polaron  method  to  whole  states  which  allowed  accurate 
band  calculations  for  insulating  solids  for  the  first  time.  This 
culminated  in  1973  in  a  theory  of  core  state  absorption  of  x-rays 
by  Insulators  Including  exclton  and  multiple  exciton  formation. 

In  the  period  1975  to  1980  he  was  a  leader  in  rigorous 
implementation  of  cluster  techniques  for  condensed  matter  systems. 
Since  1981  he  had  been  active  in  including  rigorous  quantum 


chemical  methods  to  correlate  solid  state  calculations  both  in  the 


cluster  or  the  infinite  solid  mode.  These  have  produced  useful 
and  sometimes  startling  results  for  diverse  systems  such  as  solid 
CH4  and  pure  or  defective  metals.  Currently,  he  is  fostering 
research  collaboration  between  MTU,  the  University  of  Manitoba  and 
Harwell  AERE,  England,  on  self-consistent  methods  for  studies  of 
charged  defects  in  non-metallic  solids. 
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A.  Barry  Kunz 
Administrative  Experience 

In  addition  to  remaining  active  in  research  and  teaching, 
beginning  in  1980  Dr.  Kunz  has  served  as  a  consultant  to  the  Dean 
of  Sciences  and  Arts  of  Michigan  Technological  University, 
Houghton,  Michigan.  In  1983,  1984,  1985  and  1986  he  organized  a 
Great  Lakes  Summber  Research  Workshop  at  MTU  in  Condensed  Matter 
Physics.  In  1984  a  panel  was  organized  and  chaired  for  the  U.S. 
Department  of  Energy  on  Theory  and  Simulation  of  Materials 
Structures.  During  his  period  as  Head  of  the  Physics  department, 
he  has  had  the  annual  departmental  research  funding  rise  by  a 
factor  of  10  (from  $80, 000/year  to  $1 , 000 , 000/year ) .  The  faculty 
has  increased  from  16  to  23,  the  graduate  population  from  3  to  30, 
and  a  Ph.D.  program  in  Applied  Physics  is  being  implemented.  He 
is  active  in  fostering  interdisciplinary  programs  and  has  founded 
along  with  the  departments  of  Mathematics  and  Computer  Science, 
the  Center  for  Experimental  Computation  which  is  heavily  involved 
in  research  into  parallel  computing.  He  is  also  currently  active 
in  fostering  interdisciplinary  research  into  Materials  by  Design. 
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Muhlenberg  College,  Allentown,  1962. 
Lehigh  University/  Bethlehem,  1964. 
Lehigh  University,  Bethlehem,  1966. 


Consultant,  Los  Alamos  National  Laboratory 
Consultant,  Energetics,  Inc. 

Presidential  Professor,  Head,  Physics,  Michigan 
Technological  University. 

Professor,  Physics,  University  of  Illinois  Urbana 
Champaign . 

Adjunct  Professor,  Physics,  Michigan  Technological 
University . 

Director,  Institute  for  Condensed  Matter  Studies, 
Michigan  Technological  University. 

Associate  Professor,  University  of  Illinois, 

Urbana . 

Consultant,  Photo  Products,  E.  J.  Dupont. 

Assistant  Professor,  University  of  Illinois, 

Urbana . 

Research  Analyst,  U.S.  Air  Force  Aerospace 
Research  Laboratory,  WPAFB,  Dayton. 

Research  Assistant  Professor,  Physics  and 
Materials  University  of  Illinois,  Urbana. 


1966  -  1969:  Research  Associate  and  Instructor,  Lehigh 

University,  Bethlehem. 


Honors 

1987:  Editorial  Board  Computers  In  Physics 

1987:  Michigan  Universities  Board  of  Governors  Outstanding 

Faculty  award 

1986:  MTU  Distinguished  Research  Award 

1985:  Editorial  Board  International  Journal  of  Quantum 

Chemistry 

1984:  Chairman,  U.S.  Department  of  Energy,  Panel  on  Theory  and 

Simulation  of  Defects 

Elected  Fellow  of  the  American  Physical  Society 


1976: 


Computer  Modeling 
Catalysis,  Bonded  Catalysts 

Surface  Science,  Surface  Defects,  Chemisorption 
Ionic  Solids,  Molecular  Solids,  Energetic  Solids 
Defects  in  Metals 

A.B.  KUNZ  -  Publications,  1980  to  present. 

119.  D.  C.  Reynolds,  R.  J.  Almassy,  C.  W.  Litton,  G.  Koos,  A.  B. 
Kunz  and  T.  C.  Collins.  Donor-Acceptor  Recombination 
Spectra  in  CuCl .  Phys .  Rev.  Lett.  4_4  ,  204  (1980). 

120.  G.  C.  Wepfer,  G.  T.  Surratt,  R.  S.  Weidman  and  A.  B.  Kunz. 
Theoretical  Study  of  H  Chemisorption  on  NiO,  II.  Surface 
and  Second-Layer  Defects.  Phys.  Rev.  B  2_1,  2596  (1980). 

121.  A.  B.  Kunz  and  D.  M.  Ginsberg.  Band  Calculations  of  the 
Effect  of  Magnetic  Impurity  Atoms  on  the  Properties  of 
Superconductors.  Phys.  Rev.  B  22,  3165-3172  (1980). 

122.  A.  K.  Ray,  S.  B.  Trickey,  R.  S.  Weidman  and  A.  B.  Kunz. 
Lattice  Constant  at  the  Insulator-Metal  Transition  of 
Crystalline  Xenon.  Phys.  Rev.  Lett.  ±5,  933-935  (1980). 

123.  A.  B.  Kunz,  R.  S.  Weidman,  J.  Boettger,  and  G.  Cochran. 
Theory  of  Electronic  States  in  Narrow-Band  Materials; 
Applications  to  Solid  Ar,  CuCl  and  CdS .  Int .  J.  Quant. 

Chem.  S  14,  585-596  (1980). 

124.  H.  Itoh,  G.  Ertly  and  A.  B.  Kunz.  Molecular  Orbital  Study 
on  the  Interaction  of  Single-Transition  Metal  Atoms  with  NH3 
and  H2 0  Ligands,  Z.  Naturforsch.  36  a,  347-353  (1981). 


125.  A.  B.  Kunz  and  T.  0.  Woodruff.  Some  Highly  Excited  States 
of  the  LiH  Molecule:  Calculations  Relevant  to  Core-Hole 
Initiated  Relaxation  Processes  (CHIRP's)  in  Ionic  Crystals, 
Sol.  St.  Comm.  38,  629-631  (1981). 

126.  K.  L.  Bedford  and  A.  B.  Kunz,  Ab  Initio  Studies  of  the 
Electronic  Properties  of  the  Silica  Surface,  Sol.  St.  Comm. 
38,  411-414  (1981). 

127.  A.  B.  Kunz  and  J.  T.  Waber.  Positron  and  Electron  Energy 
Bands  in  Several  Ionic  Crystals  Using  Restricted  Hartree- 
Fock  Method.  Solid  ST.  Commun.  39,  831-834  (1981). 


128.  A.  B.  Kunz,  R.  S.  Weidman  and  T.  C.  Collins.  Pressure 
Induced  Modifications  of  the  Energy  Band  Structure  of 
Crystalline  CdS.  J.  Phys.  C  14,  L581-L584  (1981). 
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129.  H.  Itoh,  G.  Ertl  and  A.  B.  Kunz .  A  Molecular  Orbital  Study 
on  the  Interaction  of  Dinitrogen  with  Transition  Metal 
Atoms.  Chemical  Physics  59,  149-156  (1981). 

130.  A.  B.  Kunz,  Electronic  Structure  of  NiO.  J.  Phys .  C  1 4 , 
L455-L460 ,  (1981). 

131.  R.  S.  Weidman,  K.  L.  Bedford  and  A.  B.  Kunz.  Ab  Initio 
Calculations  of  the  Electronic  Properties  of  Polyethylene, 
Sol.  St.  Commun.  39,  917-919  (1981). 

132.  A.  B.  Kunz,  Properties  of  Narrow  Band  Insulators,  Int .  J. 
Quant.  Chem.  515,  487-498  (1981). 

133.  J.  M.  Blaisdell  and  A.  B.  Kunz.  Theoretical  Study  of  0 
Chemisorption  on  NiO.  Solid  St.  Commun.  40,  745-747  (1981). 

134.  A.  K.  Ray,  S.  B.  Trickey  and  A.  B.  Kunz.  BDD  -  FCC 
Allotropy,  F-Bands  and  Metallization  in  Xe  and  Kr .  Solid 
ST.  Commun.  41. ,  351-353  (1982). 

135.  K.  L.  Bedford  and  A.  B.  Kunz.  Ab  Initio  Studies  of  the 
Initial  Adsorption  of  Oxygen  onto  the  Aluminum  (100) 

Surface.  Phys.  Rev.  B  25.  2119-2123  (1982). 

136.  L.  A.  Grunes,  R.  D.  Leapman,  C.  N.  Wilker,  R.  Hoffmann,  and 
A.  B.  Kunz.  Oxygen  K  near-edge  Fine  Structure:  An 
Electron-energy  Loss  Investigation  with  Comparison  to  New 
Theory  for  Selected  3d  Transition-metal  Oxides.  Phys.  Rev. 

B  25,  7157-7173  ( 1982) . 

137.  A.  B.  Kunz.  Study  of  the  Electronic  Structure  of  Twelve 
Alkali-Halide  Crystals.  Phys.  Rev.  B  26,  2056-2069  (1982). 

138.  A.  B.  Kunz.  Electronic  Structure  of  AgF,  AgCl  and  AgBr . 
Phys.  Rev.  B  r_j,  2070-2075  (  1982). 

139.  A.  D.  Zdetsis  and  A.  B.  Kunz,  Cluster  Adsorption  of  Argon  on 
Alumnina.  Phys.  Rev  B  26,  4756-4758  (1982). 


140.  A.  B.  Kunz  and  J.  T.  Waber,  Concerning  the  Trapping  of 

Positrons  in  Ionic  Solids,  in  Positron  Annihilation,  P.  <3. 
Coleman,  S.  C.  Sharma  and  L.  M.  Diana,  Eds.,  682-683  (North 
Holland,  1982). 


141.  A.  Y.  S.  Kung,  A.  B.  Kunz  and  J.  M.  Vail,  Unrestricted 
Hartree-Fock  Cluster  Analysis  of  F  and  F  Centers  in  some 
Alkali  Halides.  Phys.  Rev.  B  26,  3352-3360  (1982). 

142.  C.  E.  Brion,  I.  E.  McCarthy,  I.  H.  Suzuki,  E.  Welgold,  G.  R. 
J.  Williams,  K.  L.  Bedford,  A.  B.  Kunz  and  R.  S.  Weidman, 
Electron  Momentum  Distributions  and  Binding  Energies  for  the 
Valence  Orbitals  of  Hydrogen  Bromide  and  Hydrogen  Iodide. 

J.  of  Elect.  Spect.  and  Rel.  Phenom.  27,  83-107  (1982). 
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G.  S.  Handler  and  A.  B.  Kunz ,  The  Self-Energy  of  a  Free 
Electron  Gas.  Chem .  Phys .  Lett.  460-462  (1982). 


A.  B.  Kunz  and  J.  T.  Waber,  A  Theoretical  Study  of  the 
Binding  of  Positrons  to  Gaseous  Molecules,  in  Positron 
Annihilation.  P.  G.  Coleman,  S.  C.  Sharma  and  L.  M.  Diana, 
Eds.,  113-115  ( Horth-Hol land ,  1982). 


J.  C.  Boisvert,  A.  B.  Kunz  and  T.  0.  Woodruff,  Core  Excitons 
in  Ionic  Crystals:  I.  A  Born-Haber-type  cycle  for  the 
energy  of  the  Core  Exciton  Associated  with  the  NaLII,III 
Edge  in  NaF.  J.  Physc.  C.  J_5 ,  5033-5036  (1982). 

A.  B.  Kunz,  J.  C.  Boisvert  and  T.  0.  Woodruff,  Core  Excitons 
in  Ionic  Crystals:  II.  Calculations  of  the  Na  LII.III  Edge 
in  NaF.  J.  Physc.  C.  15,  5037-5043  (1982). 

A.  B.  Kunz  and  T.  0.  Woodruff,  Core  Excitons  in  Ionic 
Crystals:  III.  Calculations  ( including  correlations)  of 

Atomic  Energies  required  in  the  Born-Haber  Cycle  for  the  Na 
LII.III  Edge  Exciton  in  NaF.  J.  Phys.  C.  1_5,  5045-5047 
(1982)  . 


J.  C.  Boettger  and  A.  B.  Kunz,  Ab  Initio  Energy  Bands  for 
CdO .  Phys.  Rev.  B  27 ,  1359-1362  (1982). 

A.  B.  Kunz  and  C.  P.  Flynn,  A  New  Configuration  Interaction 
Method  for  Excitons  and  Interband  Processes  Applied  to  LiF . 
J.  Phys.  C.  16,  1659-1663  (1983). 

A.  B.  Kunz  and  C.  P.  Flynn,  Excitonic  Effects  in  the 
Interband  Spectra  of  Metals.  Phys.  Rev.  Lett.  50,  1524-1527 
(  1983)  . 

Russian  reprint  of  #117,  updated  for  Soviet  consumption 
(  1983)  . 

A.  B.  Kunz,  Electronic  Structure  and  Optical  Properties  of 
Solid  Methane.  Phys.  Rev.  B  28,  3465-3473  (1983). 

J.  M.  Blaisdell  and  A.  B.  Kunz,  "Theoretical  Study  of  0 
Chemisorption  on  N10,  Perfect  Surfaces  and  Cation 
Vacancies".  Phys.  Rev.  B  29,  988-995  (1984). 

A.  B.  Kunz,  "A  Series  of  Tests  of  Small  and  Medium  Scale 
Computers  Commonly  Used  for  Computations  by  Solid  State 
Theorists  and  Quantum  Chemists".  Int .  J.  Quant.  Chem.  Synp. 
17,  623-627  (1983). 

C.  P.  Flynn  and  A.  B.  Kunz,  "Recent  Results  for  Excitonic 
Processes  in  the  Spectra  of  Metals  and  Alloys".  Int.  J. 

Quat .  Chem.  Synp.  1_7,  573-582  (1983). 


156.  A.  K.  Ray,  S.  B.  Trickey  and  A.  B.  Kunz ,  "Optical  Adsorption 
of  Solid  Xenon  at  High  Pressure".  Phys .  Stat.  Sol.  B  121, 
K47-K50  (1984). 

157.  D.  R.  Beck  and  A.  B.  Kunz,  "Excitation  Energies  for  the 
Lowest  Triplets  and  Singlet-triplet  Splittings  in  Gas-phase 
Methane  Including  Many-body  Effects".  J .  Phys.  B  .17,  2159- 
2168  (1984). 

158.  A.  B.  Kunz,  "Theory  of  the  Electronic  Structure  and  Optical 
Properties  of  Organic  Solids:  Collective  Effects",  in 
Quantuum  Chemistry  of  Polymers ,  Solid  State  Aspects ,  ed .  J . 
Ladik,  83-89  (D.  Reidel  1984). 

159.  A.  B.  Kunz,  "Chemisorption  and  the  Electronic  Structure  of 
Transition  Metal  Oxides  and  Transition  Metals  Bonded  to 
Oxide  Surface".  Phil.  Mag.  51,  209  (1985). 

160.  A.  B.  Kunz,  "Defect  Simulations  and  Supercomputers". 

Reports  on  Materials  Science,  submitted. 

161.  A.  B.  Kunz,  J.  C.  Boisvert  and  T.  0.  Woodruff,  "LiK-edge 

Soft  X-ray  Absorption  in  LiF" .  Phys.  Rev.  B30 ,  2158 

(1984). 

162.  M.  H.  Bakshi,  G.  A.  Denton,  C.  P.  Flynn,  J.  C.  Boisvert,  and 
A.  B.  Kunz,  "Low-energy  Excitonic  Resonances  in  Metals  I: 
Experiments  on  Divalent  Atoms".  Phys.  Rev.  B,  Vol .  31,  pp. 
4972-2983. 

163.  J.  C.  Boisvert,  P.  W.  Goalwin,  A.  B.  Kunz,  M.  H.  Bakshi,  and 
C.  P.  Flynn,  "Low-energy  Excitonic  Resources  in  Metals  II: 
Cluster  Calculation  for  Divalent  Impurities  in  Li".  Phys. 
Rev.  B.  Vol.  31,  pp.  4984-4992. 

164.  A.  B.  Kunz,  1985,  "The  Electronic  Structure  of  Impurities 
and  Defects  in  Simple  Metals".  Characterization  and 
Behavior  of  Materials  with  Submicron  Dimensions,  World 
Scientific  Publishing  Co. 

165.  A.  B.  Kunz,  "Solid  State  Theory".  The  Encyclopedia  of 
Physics,  ed.  Besancon. 

166.  A.  B.  Kunz,  1985,  "Defect  Simulation  and  Supercomputers." 
Materials  Science  Forum,  Vol.  4,  pp .  155-164. 

167.  M.  H.  Bakshi,  G.  A.  Denton,  C.  P.  Flynn,  J.  C.  Boisvert  and 
A.  B.  Kunz,  "Low-Energy  Excitonic  Resonances  in  Metals,  I. 
Experiments  on  Divalent  Atoms",  Physical  Review  B. ,  Vol.  31, 
pp.  4972-4983.  1985. 
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J.  C.  Boisvert,  P.  W.  Goalwin,  A.  B.  Kunz ,  M.  H.  Bakshi  and 
C.  P.  Flynn,  "Low-Energy  Excitonic  Resonances  in  Metals,  II. 
Cluster  Calculat ions  for  Divalent  Impurities  in  Li  Metal”, 
Physical  Review  B.,  Vol .  31,  pp .  4984-4992.  1985. 

N.  C.  Bacalis  and  A.  B.  Kunz,  1985,  "Ab  Initio  Calculations 
of  Selected  Ionization  Spates  of  Cu  on  MgO"  (001),  Physical 
Review  B,  Vol.  32  pp.  4857-4865. 

A.  D.  Zdetsis  and  A.  B.  Kunz,  1985,  "Ab  Initio  Cluster  Study 
of  Ni  Adsorption  on  Alumina",  Physical  Review  B,  Vol.  32, 
pp.  6358-6362. 

A.  B.  Kunz,  C.  P.  Flynn,  J.  C.  Boisvert,  M.  Bakshi, 
"Electronic  Structure  of  Impurtities  and  Defects  in  Sample 
Metals,  in  Characterization  and  Behavior  of  Materials  with 
Submicron  Diments ions " ,  World  Scientific  Publications,  pp. 
83-100,  1986. 

A.  B.  Kunz,  "Strategies  of  Studying  Point  Defects", 
Supercomputers ,  edited  by  F.A.  Matsen  and  T.  Tajima, 
University  of  Texas  Press,  pp.  250-268,  1986. 

A.  B.  Kunz,  "Computational  Consideration  in  Solving  the  Many 
Electron  Problem  for  Molecular  System",  Applications  of 
Supercomputers ,  edited  by  D.F.  Lockhart  and  D.L.  Hicks, 
Journal  of  Applied  Mathematics  and  Computation,  pp.  23-33, 
1986. 

P.  W.  Goalwin,  A.  B.  Kunz,  "Perturbation-Theory  Method  of 
Calculating  the  Energies  and  Excitation  Energies  of  Atomic, 
Molecular  and  Solid  State  Systems",  Physical  Rev.  B  34,  pp. 
2140-2146,  1986. 

P.  W.  Goalwin  and  A.  B.  Kunz,  "Ab  Initio  Calculations  of  the 
Excitonic  Spectra  of  the  Alkali  Halides  Including 
Relativistic  Effects",  Phys.  Rev.  B35,  5795  (1987). 

A.  B.  Kunz,  "Theoretical  Study  of  Defects  and  Chemisorption 
by  Oxide  Surfaces,  for  Theory  of  Oxide  Surfaces  and 
Properties",  J.  Nowtny  Ed.,  in  press. 

C.  P.  Flynn  and  A.  B.  Kunz,  "Comment  on  Role  of  Electron- 
Hole  Interactions  in  the  Crystal  Spectra  of  Metals",  Phys. 
Rev. ,  in  press . 

A.  B.  Kunz,  "A  Cluster  Model  for  the  Adsorption  of  Rare  Gas 
Atoms  by  Simple  Metals",  submitted. 

A.  B.  Kunz  and  J.  M.  Vail,  "Quantum  Mechanical  Cluster 
Lattice  Interaction  in  Crystal  Simulation",  submitted. 

A.  B.  Kunz  and  D.  R.  Beck,  "Possible  Role  of  Charged  Defects 
in  Molecular  Solids",  Phys.  Rev.,  in  press. 
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181.  A.  B.  Kunz,  J.  Meng  and  C.  Woodward,  "Computational 

Consideration  for  the  Theoretical  Study  of  Defects  in  Solid 
State  Systems",  submitted. 

162.  A.  B.  Kunz,  "Embedded  Cluster  Method  for  Surface  and 
Chemisorption  Studies",  submitted. 


A.  B.  KUNZ  -  Invited  Papers 


1969  29  September  to  4  October,  European  Congress  on  Molecular 

Spectroscopy,  "Energy  Band  Structure  and  Optical  Properties 
of  Four  Alkali  Fluorides",  Universite'  de  Liege,  Belgium. 

1972  4  July,  International  Conference  on  Electron  Energy 
STructure,  "Ab  Initio  Energy  Band  STudies  for  Ionic  Solids", 
University  of  Exeter,  England. 

1973  18-30  June,  NATO,  Advanced  Study  Institute,  On  Elementary 
Excitation  in  Solids,  Molecules  and  Atom.  Two  Lectures: 

"Ab  Initio  Energy  Band  Methods" 

"Soft  X-ray  Spectra  of  Molecules". 

27-29  December,  APS  winter  meeting  at  Berkeley,  "Soft  X-ray 
Spectra  of  Insulating  Solids",  University  of  Berkeley, 
California. 

1974  January,  International  Symposium  on  Quantum  Theory  of  Atoms, 
Molecules,  and  the  Solid-State.  "Many-Body  Effects  in 
Energy  Band  Theory",  Sanibel  Island,  Florida  (University  of 
Florida) . 

22-26  July,  IV  International  Conference  on  Vacuum 
Ultraviolet  Radiation  Physics,  "Calculation  of  Optical 
Absorption  Edges  and  Photo-Emission  Edges  in  Solids", 
Hamburg,  Germany. 

17-18  November,  Midwest  Solid  State  Theory  Symposium, 
"Photo-Emission  Studies  of  Insulating  Solids",  Indiana 
University. 

December,  Peat  Symposium  of  the  U.S.  Geological  Society,  in 
cooperation  with  T.  L.  Phillips  and  D.  J.  Mickish, 
"Paleobotany  of  Petrified  Peat  (Coal  Ball)  Coal  Member  of 
the  Illinois  Basin",  Miami  Beach,  Florida. 

1975  23-27  June,  Summer  Research  Conference  in  Theoretical 
Chemistry,  Session  Chairman  on  Subject  of  Ab  Inition  Theory 
of  Crystal  and  Surface  Science,*  Boulder,  Colorado. 

15  July  to  15  September,  "Centre  European  de  Calcul  Atomique 
et  Maelculaire,  Workshop  on  Ab-initio  one-electron 
potentials".  Seminar  series  declined,  Orsay,  France. 

1976  18-24  January,  International  Symposium  on  the  Quantum  Theory 
of  Atoms,  Molecules  and  the  Solid-State;  Symposium  Leader; 
Sanibel  Island,  Florida. 

9-13  August,  Gordon  Conference  on  Atomic  and  Molecular 
Interactions . 


1977  5  November,  Midwest  Solid  State  Theory  Meeting,  Argonne  Nat. 
La.,  declined  due  to  conflict. 

1978  4-18  June,  NATO  ASI,  Excited  States  in  Quantum  Chemistry, 
KOS,  Greece. 

1979  17  March,  Int.  Symp.  on  Quantum  Chemistry,  Electronic 
Properties  of  CuCl,  Palm  Coast,  Florida. 

25  August,  Sagamore  Conference  on  Crystallography,  Mt . 
Tremblant,  Canada. 

1980  15  March,  International  Sym .  on  Quantum  Chemistry,  Recent 
Developments  in  CuCl  and  CdS,  Palm  Coast,  Florida. 

22  March,  IBM  Workshop  on  Effective  One-Electron  Potentials. 
Electronic  Structure  of  Transition  Metal  Oxides,  Ossining, 
New  York. 

1-4  July,  Gordon  Conference  on  Electron  Distribution  and 
Chemical  Bonding  -  discussion  leader. 

1981  21  January,  Workshop  on  Fundamental  Research  on  Energetic 
Materials,  Finite  Cluster  Techniques  for  Modeling  Extended 
Systems,  Berkeley,  California. 

14  March,  International  Symposium  on  Quantum  Chemistry,  Palm 
Coast,  Florida.  Fundamental  Properties  of  Narrow  Band 
Materials . 

June,  Symposium  of  NRCC,  Los  Alamos.  Two  invited  talks: 

1.  Simulation  of  Extended  Systems  by  Finite  Clusters 

2.  Solid  State  Effective  Potentials 

1982  May,  DOD  Workshop  on  Energetic  Materials:  Properties  of 
Energetic  Solids  and  Molecular  Crystals. 

1983  12  March,  Sanibel  Symposium  on  Impact  of  Computers  on 
Quantum  Chemistry  and  Solid  State  Computations  -  Panel. 

24  July  -  5  August,  "Excited  States  in  Molecular  Solids, 
Collective  Effects",  NATO  ASI,  Braunloge,  West  Germany  on 
Quantum  Theory  of  Polymers,  Solid  State  Aspects. 

9  August,  "Energy  Transfer,  Trapping  and  Release  in 
Molecular  Crystals  and  Energetic  Solids",  DOD  Workshop  on 
Energetic  Materials,  Great  Oaks,  Maryland. 

3  October,  "The  Electronic  Structure  of  Impurities  and 
Defects  in  Simple  Metals",  Fall  Meeting  of  the  Metallurgical 
Society  and  the  American  Society  of  Metals. 

1984  1-10  March,  Sanibel  Symposium,  Parts  I  and  II,  Organizing 
Committee  for  Part  I  and  Session  Chairman. 


9-13  April,  Royal  Society  of  Chemistry,  London,  England; 
Invited  Plenary  Lecture:  "Chemisorption  and  the  Electronic 
Structure  of  Transition  Metal  Oxides  and  Transition  Metals 
Bonded  to  Oxide  Surfaces". 

6-10  August,  Department  of  Energy  Panel  on  Simulation  and 

Theory  of  Defect  Structures  in  Solids  - ,  Chairman  and 

organizer . 

5-8  September,  Midwest  Solid  State  Theory  Conference, 
Minneapolis,  Minnesota,  "Defect  Simulation  and  Super 
Computers" . 

1985  17-20  March,  High  Speed  Computers  and  the  Future  of 
Scientific  Computation,  University  of  Texas,  Austin,  TX 
"Theory  of  Defects  and  Supercomputers". 

24-29  June,  DoD  workshop  on  Initiation  Fundamentals,  Great 
Oak,  MD,  "Defects  in  Energetic  Solids". 

5-6  August,  Super  Computer  Conference,  MTU,  Houghton,  MI, 
"Strategies  For  Studying  Point  Defects". 

1986  29-31  January,  DoD  Workshop  on  Energetic  Solids,  Silver 
Spring,  MD  NSWC,  "Defects  and  Initiation". 

12-16  March,  Sanibel  Symposium,  Marineland,  FL,  Session 
Chairman . 

February  1986,  D.O.E.  Workshop  on  Energetics,  White  Oak, 
Maryland,  "Mechanism  for  Destabilization  of  Engergetic 
Solids". 

October  1986,  D.O.E.  Conference  on  Initiation,  Los  Alamos, 
New  Mexico. 

1987  February  1987,  FPS  Symposium  on  Parallel  Computer,  "Computer 
Modeling  of  Solid  State  Defects",  Dallas,  Texas. 

March  1987,  Sanibel  Symposium,  "Computational  Consideration 
for  Solving  the  Many  Body  Problem  in  Solids",  Marineland, 
Florida. 

May  1987,  DoD  Workshop,  "Theory  of  Defects  in  Energetic 
Solids",  Great  Oak,  Maryland. 

March  1987,  APS  March  Meeting  Pre  Symposium,  invitation 
declined  due  to  conflict. 

March  1987,  DoD-DoE  Symposium  on  Propellants,  invitation 
declined  due  to  conflict. 


August  1987,  Wright-Patterson  Symposium  on  Advanced 
Computing,  "Theory  of  Defects  in  Solids:  Computational 
Strategies,  Dayton,  Ohio. 

‘These  sessions  consist  of  a  few  short  invited  talks  and  an 

extemporized  program. 

Seminars  given  at  the  following: 

National  Bureau  of  Standards 

Aerospace  Research  Laboratory,  Wright-Patterson  AFB 

Brookhaven  National  Laboratory 

Argonne  National  Laboratory 

Oak  Ridge  National  Laboratory 

Bell  Telephone  Laboratories  (Murray  Hill) 

Bell  Telephone  Laboratories  (Holmdale) 

Eastman-Kodak  Comapny 

E.  I.  Dupont  de  Nemours  Company 

Belgian  Atomic  Energy  Establishment-Moll 

University  of  Liege 

Lehigh  University 

Drexel  Institute  of  Technology 

Battelle  Memorial  Institute 

Queens  College  CUNY 

University  of  North  Carolina 

Wake-Forrest  University 

University  of  Georgia 

Stanford  University 

Arizona  State  University 

Oklahoma  State  University 

University  of  Illinois 

Indiana  University 

Purdue  University 

Louisiana  State  University  (Baton  Rouge) 

University  of  Manitoba 
Texas  A  &  M  University 
Naval  Research  Laboratories 

U.  S.  Army  Ballistics  Research  Laboratory,  Aberdeen,  Maryland 

Michigan  State  University 

Hirman  College 

General  Motors  Corporation 

Michigan  Technological  University 

Air  Force  Avionics  Laboratory 

University  of  Missouri,  Columbia 

University  of  Missouri,  St.  Louis 

University  of  Missouri,  Rolla 

Los  Alamos  National  Laboratory 

University  of  New  Orleans 

Johns  Hopkins  University 

University  of  Florida 

Western  Kentucky  University 

Washington  State  University 

Western  Kentucky  University 


VSVWWV  Vvwvm.  ¥V"W-Trtriru 


-  JTj r»jr«jr» it- 


( seminar ) 
NLTNIF 


Washington  State  University 
Exxon  Corporation,  Clinton,  NJ 
Oak  Ridge  National  Laboratory, 

Spokane,  WA,  APS  Conference 

San  Diego,  CA,  DARPA  Conference  on  Materials  Science 
Los  Alamos,  NM,  LANL-BES  Review 
Albuquerque,  NM,  Sandia,  BES  Review 

Boston,  MA,  Winter  Mtng,  Materials  Research  Society, 
Session  Chair,  Editor 
State  University  of  New  York,  Albany 
U.  S.  Naval  Research  Laboratory 
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1971-1976 


1976-1981 


1974-1985 


A.  Barry  Kunz 

Funding  History  (Major  Grants) 


Contract  for  "Theory  of  Correlation  in  Solids", 
Aerospace  Research  Laboratory,  Wright-Patterson  AFB 
Initial  annual  level  ~$  10,000 

Final  annual  level  ~$  50,000 

Total  award  "$150,000 


Contract  for  "Band  Theory  of  Semiconductiors" ,  Air 
Force  Office  of  Scientific  Research. 

Initial  annual  level  ~$  50,000 
Final  annual  level  ~$  70,000 
Total  award  "$300,000 


Grant  for  "Theory  of  Solids",  National  Science 
Foundation  via  University  of  Illinois  Materials 
Research  Laboratory. 

Initial  annual  level  ~$  70,000 

Final  annual  level  "$180,000 

Total  award  "$1,210,000 


Grant  for  "Energy  Trapping,  Release,  and  Transport  in 
Energetic  Solids",  Office  of  Naval  Research. 

Initial  annual  level  ~$  50,000 
Final  annual  level  "$  80,000 
Total  award  $346,000 


1984-Present  Grant  for  "Theory  of  Defects",  Department  of  Energy 

Current  annual  level  $  96,000 

Total  project  to  date  $362,000 


1985- Present  Grant  for  "Defects  in  Non  Metals",  Office  of  Naval 

Research . 

Current  annual  level  $278,000 

Total  project  to  date  $816,000 

1986- Present  Contract  for  "Science  Support  for  Process  Rheology", 

Office  of  Naval  Research,  subcontract  from  Lockheed, 
Inc . 

Current  annual  level  $150,000 

Total  project  to  date  $300,000 


1981-1986 


1 


VITA 


Bryan  H.  Suits 

Bryan  H.  Suits  was  born  in  Ann  Arbor,  Michigan  in  the  year 
1954.  While  in  High  School  he  received  several  awards  for 
academic  achievement  in  math  and  physics  including  a  Navy  Science 
Award,  an  Exxon  "Scientist  of  Tomorrow"  award  as  well  as  an  award 
from  the  state  of  Michigan  for  his  performance  on  the  state  wide 
mathematics  contest.  He  graduated  from  Ann  Arbor  Huron  High  school 
with  honors  in  1973. 

Dr.  Suits  continued  his  education  at  Carleton  College  in 
Northfield,  Minnesota.  During  this  time  he  spent  a  six  month 
internship  at  Argonne  National  Labs  studying  rare  earth  hydrides 
using  Mossbauer  spectroscopy  with  Gopal  Shenoy  and  Bobby  Dunlap. 
Upon  returning  to  Carleton  he  designed  and  completed  a  senior 
project  studying  gamma  ray  angular  correlations  in  Pd.  He 
obtained  his  BA  with  a  major  in  physics  in  1977  magna  cum  laude 
and  was  elected  to  the  honorary  fraternities  Phi  Beta  Kappa  and 
Sigma  Xi . 

From  the  spring  to  fall  of  1977  he  worked  at  the  Environmental 
Research  Institute  of  Michigan  (ERIM)  in  Ann  Arbor.  This  work 
centered  on  the  use  of  computer  modeling  to  establish  criteria  for 
the  aerial  evaluation  and  identification  of  vegetation. 

In  the  fall  of  1977  he  entered  graduate  school  at  the 
University  of  Illinois  at  Urbana-Champaign .  In  the  summer  of  1978 
he  began  work  with  Prof.  C.  P.  Slichter  studying  charge  density 
waves  in  transition  metal  dichalcogenides  using  NMR .  This  work 
included  growing  and  characterizing  crystalline  samples  as  well  as 
NMR  and  NQR  measurements.  In  January  of  1979  he  was  granted  a 
Master's  of  Physics.  His  Ph . D  thesis,  titled  "Confirmation  of 
McMillan's  Concept  of  Discommensurat ions  in  2H-Tantalum 
Diselenide"  was  submitted  and  accepted  in  August  of  1981.  From 
September  1981  to  the  fall  of  1982  he  worked  as  a  research 
associate  with  Slichter 's  group.  During  this  time  he  designed  and 
performed  NQR  experiments  on  the  novel  conductor  niobium 
triselenide . 

In  the  late  fall  of  1982  Dr.  Suits  began  work  at  the 
University  of  Pennsylvania  as  a  research  associate  with  the 
Laboratory  for  Research  on  the  Structure  of  Matter  with  Profs.  G. 
Farrington  and  David  White.  This  work  emphasized  the  study  of 
ionic  motion  in  the  0-alumina  family  of  materials  using  NMR. 

While  at  Pennsylvania,  Dr.  Suits  began  his  pioneering  work  on  the 
use  of  NMR  imaging  to  study  solids.  In  particular,  producing 
Images  of  diffusing  ions  in  the  "superionic  conducting"  0-alumina 
materials.  He  also  developed  and  taught  an  experimental  cross- 
disiplinary  course,  "Digital  Circuitry  and  Design  of 
Instrumentation",  a  graduate  level  special  topics  course.  This 
course  was  the  first  of  its  kind  to  be  sponsored  by  the  LRSM. 


In  the  fall  of  1985,  Dr.  Suits  began  his  current  appointment 
as  an  assistant  professor  with  the  department  of  physics  at 
Michigan  Technological  University.  He  is  developing  non-invasive 
HMR  imaging  techniques  to  study  damaged  solids  and  to  study 
complex  fluid  flow  in  particle  filled  polymer  systems.  This  work 
is  being  supported  by  two  contracts  with  the  Office  of  Naval 
Research,  the  second  of  which  is  a  joint  contract  with  Lockheed 
Corporation.  Dr.  Suits  is  a  member  of  the  American  Physical 
Society  and  the  Society  for  Magnetic  Resonance  Imaging. 

Publications 

"Valance  Transformation  of  Eu  Atoms  Due  to  Defect  Association  in 
Cold  Worked  EuMg  Alloys,"  R.  G.  Pi  rich,  C.  R.  Burr,  G.  K. 
Shenoy,  B.  D.  Dunlap,  B.  Suits  and  J.  D.  Phillips,  Phys .  Rev . 
Lett.  38,  1142  ( 19) .77 

3  + 

"Mossbauer  Study  of  Er  in  ErH_ , "  B.  Suits,  G.  K.  Shenoy,  B.  D. 
Dxinlap,  and  D.  G.  Westlake,  J.  of  Magnetism  and  Magnetic 
Materials  5,  344  (1977). 

"Magnetic  Properties  of  Dihydrides  and  Dideuterides  of  Er  and  Ho," 
J.  M.  Friedt,  B.  Suits,  G.  K.  Shenoy,  B.  D.  Dunlap,  and  D.  G. 
Westlake,  J.  Appl .  Phys.  50,  2049  (1979). 

"Observation  of  Newly  Discovered  112K  Charge  Density  Wave 

Transition  in  2H-TaSe  Using  Electrical  Resistance”,  B.  H. 

Suits  and  M.  C.  Chen,  Phys .  Let .  79A .  224  (1980). 


"Confirmation  of  McMillan's  Concept  of  Discommensurations",  B.  H. 
Suits,  S.  Couturie,  and  C.  P.  Slichter,  Phys .  Rev.  Lett .  45, 
194  (1980). 

"NMR  Test  of  McMillan's  Concept  of  Discommensurations  in 

2H-TaSe_,"  B.  H.  Suits,  S.  Couturie,  and  C.  P.  Slichter,  Phys . 
Rev.  B23,  5142  (1981) . 

"Thermopower  of  TaS",  R.  L.  Allgeyer,  B.  H.  Suits,  and  F.  C. 
Brown,  Solid  State  Commun .  43,  207  (1982). 

"An  NQR  Study  of  the  145K  CDW  Transition  in  NbSe-" ,  B.  H.  Suits 
and  C.  P.  Slichter,  Phys.  Rev.  B29,  41  (1984)7 

"NMR  Imaging  in  Solids",  B.  H.  Suits  and  David  White,  Solid 
State  Commun .  50,  291  (1984). 

"Chemical  Shifts  of  Lead  Halides",  B.  H.  Suits  and  David  White, 

In  preparation. 

"Nuclear  Magnetic  Resonance  Imaging  of  Temperature  Profiles",  B.H. 
Suits  and  David  White,  Appl .  Phys .  60,  3772  (1986). 
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"A  Nuclear  Magnetic  Resonance  Study  of  Al-Mn  Quasicrystals  and 
Related  Materials",  K.  R.  Carduner,  B.  H.  Suits,  J.  A.  DiVerdi ,  M. 
p.  Murphy  and  David  White,  J.  Mater.  Res.  2,  431  (1987). 

"Imaging  Electromagnetic  Fields  Using  the  Magnetic  Resonance  of  a 
Paramagnetic  Gas",  S.  A.  Marshall,  B.  H.  Suits,  M.  T.  Umlor,  and 
U.  N.  Zhang,  to  be  published  in  J.  Magn.  Reson. 

"An  NMR  Study  of  Li  in  AlLiCu  Icosahedral  Alloys",  C.  Lee,  David 
White,  B.  H.  Suits,  P.  A.  Bancel,  and  P.  A.  Heiney,  submitted  to 
the  Physical  Review. 

Reports 

J.  Robert  Maxwell  and  Bryan  Suits,  ERIM  reports  127600-7-f  (190 
pgs)  and  127600-3-f  (134  pgs )  for  the  U.S.  Dept,  of  Justice, 
Contract  DEA-77-6. 

Significant  Presentations 

"NMR  Confirmation  of  McMillan's  Concept  of  Discommensurat ions" , 
Invited  talk  for  the  Division  of  Condensed  Matter  Symposium  on 
CDW's  and  SDW's  in  Incommensurate  Systems,  March  Meeting  of  the 
American  Physical  Society,  Phoenix  (1981). 

"NMR  Imaging  in  Solids",  Invited  talk  for  the  FACSS  fall  meeting 
in  Philadelphia  (1985). 
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Biographical  Sketch. 

S.  A.  Marshall,  Jr. 

Present 
Position : 

Telephone : 

Mailing  Address: 


Education : 


Professional 
Experience  : 


Research 
Interests : 


Organization 

Membership: 


Additional 

Professional 

Activities 


Professor,  Department  of  Physics 
Michigan  Technological  University. 

906/487-2900  (work),  906/482-5242  (home) 

Department  of  Physics,  Fisher  Hall 
Michigan  Technological  University 
Houghton,  Michigan  49931. 

Ph.D.  (physics)  The  Catholic  University  of 
America,  Washington,  D.C.,  M.S.  (physics) 
The  University  of  Michigan,  Ann  Arbor, 
Michigan,  B.S.  (physics)  The  Illinois 
Institute  of  Technology,  Chicago, 

Illinois. 


1981-  ,  Professor  Department  of  Physics, 
Michigan  Technological  Unverisity, 
Houghton,  Michigan. 

1963-1981,  Research  Scientist,  Solid  State 
Sciences  Division,  Argonne  National 
Laboratory,  Argonne,  Illinois. 

1956-1963,  Senior  Research  Scientist,  IIT 
Research  Institute,  Chicago,  Illinois. 
1952-1956,  Physicist,  U.S.  Naval  Ordnance 
Laboratory,  White-Oak,  Silver  Spring,  MD. 
1951-1952,  Physicist,  National  Bureau  of 
Standards,  Washington,  D.C. 

Physics  of  Solids:  Magnetic  Defects  in 
Nearly  Perfect  Crystals;  Paramagnetic 
Resonance  Absorption  Spectroscopy  of 
Transition  Metal  Ions  and  Free-Radical 
Molecule  Ions  in  Solids  and  Liquids; 
Nuclear  Magnetic  Resonance  Spectroscopy  in 
Hydrides  and  Superionic  Conductors: 

Crystal  Field  Theory;  Hyperfine 
Interactions;  Color  Centers  in  Alkali 
Halide  Crystals. 

The  American  Physical  Society  (1952), 

The  Society  of  Sigma  Xi  (1956), 

A.A.A.S.  (1956), 

The  Physics  Club  of  Chicago  (1958). 

1956-57,  Visiting  Scientist  at  the  U.S. 
Naval  Ordnance  Laboratory  where  some 
design  work  was  done  on  a  high  temperature 


I 


microwave  spectrometer.  1957-62,  On 
various  occasions  during  this  period 
taught  graduate  and  undergraduate  courses 
at  the  Illinois  Institute  of  Tech.  1960- 
63,  Directed  one  M.S.  and  two  Ph.D. 
thesis  programs  at  the  Ill.  Inst,  of  Tech. 

1964,  Attended  a  summer  institute  at  the 
University  of  Colorado  (Boulder)  on  the 
Theory  of  Groups  and  on  the  Lorentz  Group. 
1965-present,  Taught  a  course  on  Electron 
Spin  Resonance  Absorption  Spectroscopy  at 
a  NATO  summer  school  in  Freiburgh,  W. 
Germany;  directed  the  research  of 
approximately  20  honors  students  at  the 
average  of  two  students  per  year;  referee 
of  papers  to  scientific  journals  such  as 
The  Physical  Review,  The  Journal  of 
Magnetic  Resonance,  The  Journal  of  Physics 
and  Chemistry  of  Solids,  The  Journal  of 
Physical  Chemistry,  and  others  on 
occasion . 

For  about  seven  years,  held  position  of 
adjunct  professor  of  physics  at  the 
Northern  Illinois  University  where  an  M.S. 
thesis  was  directed  during  the  last  two 
years  of  the  appointment. 

During  the  1970-71  participated  in  a 
program  of  research  with  Dr.  S.  V.  Nistor 
of  the  Institute  of  Atonic  Physics, 
Bucharest,  Romania.  The  program  was 
sponsored  in  part  by  the  National  Science 
Foundation  (U.S.A.),  the  Ministry  of 
Science  and  Culture  (Romania)  and  the  U.S. 
Atomic  Energy  Commission. 

Academic  1960,  Irving  Segal,  M.S.  in  physics  at  the 

Research  Illinois  Institute  of  Technology. 

Supervision:  1962,  Dickron  Mergerlan,  Ph.D.  in  physics 

at  the  Illinois  Institute  of  Technology. 

1965,  Raymond  R.  Serway,  Ph.D.  in  physics 
at  the  Illinois  Institute  of  Technology. 
1970,  Joachim  Doehler,  M.S.  in  physics  at 
the  Illinois  Institute  of  Technology. 


Publications : 


36.  Temperature  Dependence  of  the  Hyperfine  Structure  Splitting 
of  Divalent  Manganese  in  Single  Crystal  Calcite. 

C.  Y.  Huang,  R.  S.  Kent  and  S.  A.  Marshall,  Phys .  Rev.  B7 , 
552  (  1973)  . 


37.  Temperature  Dependence  of  the  Hyperfine  Structure  Splittings 
of  Some  Trioxide  Molecule  Ions  in  Single  Crystal  Calcite. 
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